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Abstract 

These lectures introduce some of the basic methods of perturbative QCD and their 
, apphcations to phenomenology at high energy. Emphasis is given to techniques that are 

used to study QCD and related field theories to all orders in perturbation theory, with 
' introductions to infrared safety, factorization and evolution in high energy hard scattering. 
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1 Introduction 



o 

:a 

Quantum chromodynamics is the component of the Standard Model that describes the strong 
^ ■ interactions. It is a theory formulated in terms of quarks and gluons at the Lagrangian level, 
but observed in terms of nucleons and mesons in nature. Partly because of this, its verification 
^ ! required formulating new methods and asking new questions on how to confront quantum field 
^ [ theories with experiment, especially those that exhibit strong coupling. In what might seem an 
almost paradoxical outcome, new roles were recognized for perturbative methods based on the 
elementary constituents of the theory. 

Since strong coupling is a feature of many extensions of the Standard Model and of string 
theories, these theoretical developments in perturbative QCD remain of interest in a more general 
context. In addition, the perturbative analysis of hadronic scattering and final states played an 
essential role in the verification of the electroweak sector of the Standard Model, through the 
discoveries of the W and Z bosons and of the bottom and top quarks. Likewise, perturbative 
methods are indispensable in the search for new physics at all colliders, which will keep them 
relevant to phenomenology for the foreseeable future. Finally, at lower but still substantial 
energies, new capabilities to study the polarization-dependence of hadronic scattering and the 
strong interactions at high temperatures and densities are opening new avenues for the study of 
quantum field theory. 

^To appear in proceedings of the 2004 Theoretical Advanced Study Institute in Elementary Particle Physics 
(TASI 2004), University of Colorado at Boulder, June 6 - July 2, 2004. 
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These lectures, of course, cannot cover the full range of QCD, even at high energies in per- 
turbation theory. They are meant to be an introduction for students familiar with the basics of 
quantum field theory and the Standard Model, but not with applications of perturbative QCD. 

1 will emphasize the formulation and calculation of observables using perturbative methods, 
and discuss the justification of perturbative calculations in a theory like QCD, with csssentially 
nonperturbative long-distance behavior. 1 will also introduce elements of the phenomenology 
of perturbative QCD at high energies, keeping in mind its use in the search for signals of new 
physics. 

2 QCD from Asymptotic Freedom to Infrared Safety 

We begin with a quick prehistory of quantum chromodynamics, and of the strands of experiment 
and theory that converged in the concept of asymptotic freedom. We also show how infrared 
safety makes possible a wide range of phenomenological implementations of asymptotic freedom. 

2.1 Why QCD? 

The search that culminated in the SU (3) gauge theory quantum chromodynamics was triggered 
by the discovery of an ever-growing list of hadronic states, including mesons, excited nucleons 
and a variety of resonances, beginning with the A resonances first seen in the early 1950s [1], 
and continuing to this day [2]. Already in 1949, the idea was put forward that mesons might 
be composites of more elementary states (of nucleons in the very early work of Fermi and Yang 
[3]). Clearly hadrons were not themselves elementary in the normal sense, but for a long time 
the strong couplings between hadronic states, and their ever-growing number, left the situation 
unsettled. 

In broad terms, two alternative viewpoints were entertained for nucleon and other hadronic 
states. In bootstrap scenarios, each hadronic state was at once deducible from the others, yet 
irreducible, in the sense that no hadronic state was more elementary than any other. This idea 
eventually evolved through the dual model [4] into string theory, leaving the strong interactions 
behind, at least temporarily. The more conservative viewpoint regarded hadrons as composite 
states of elementary particles, but the observed couplings seemed to pose an insurmountable 
barrier to a field-theoretic treatment. Still, the growing success of quantum electrodynamics held 
out at least the possibility that such a picture might somehow emerge. 

With the weak and electromagnetic interactions serving as diagnostics, what we now call 
fiavor symmetry was developed out of the spectroscopy of hadronic states. Their interactions 
through local currents also hinted at point- like substructure. The weak and electromagnetic 
currents, in their turn, were naturally described as coupling to fundamental objects, and the 
quark model was born, partly as a wonderfully compact classification scheme, and partly as a 
dynamical model of hadrons. 

Increases in accelerator energies eventually clarified the situation. A prologue, both scien- 
tifically and technologically, were the 1950's measurements of electron-proton elastic scattering 
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[5] . Imposing the relevant symmetries of the electromagnetic and strong interactions (parity and 
time-reversal invariance), the cross section for electron- nucleon unpolarized scattering at fixed 
angle is customarily expressed in the nucleon rest frame, as 
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The term in square brackets is the (Mott) cross section for the scattering of a relativistic electron 
of energy Ehj a fixed source of unit charge, while E' is the outgoing electron energy. The variable 
T = Q'^/AMjf, with Mn the nucleon mass and Q'^ — —EE'sm^{9/2) the invariant momentum 
transfer. The functions Ge{Q) and Gm{Q) are the electric and magnetic form factors. They fall 
off rapidly with momentum transfer, roughly as 

where //q ~ 0.7GeV. This simple behavior is anaolgous to the form factors of nonrelativistic 
bound states with extended charge distributions [6] . This suggests a substructure for the nucleon, 
but leaves its nature mysterious. Indeed, much is still being learned about nucleon form factors 
[7] and their relation to the quarks and gluons of QCD [8, 9]. The experiments arc difficult, 
however, because of the form factors' rapid decrease with Q. This is simply to say that when a 
proton is hit hard, it is likely to radiate or to break apart. 

As available electron energies increased it became possible by the late 1960's to go far beyond 
elastic scattering, and to study highly ("deep") inelastic reactions. In these experiments, nature's 
choice became manifest. 

The idea is simple, and analogous in many ways to the famous Rutherford experiments that 
revealed the atomic nucleus. A high energy electron beam is made to collide with a nucleon 
target. One observes the recoil of the electron, measuring the total cross section as a function of 
transfered energy and scattering angle, as indicated in Fig. 1. Although each individual event is 
generally quite complicated, a simpler pattern emerges when events at a given momentum transfer 
are grouped together. Imposing again the relevant symmetries, a general deep-inelastic scattering 
(DIS) cross section can be written in terms of dimensionless structure functions, Fi{x,Q'^) and 
F2{x, (5^), where the variable x is defined as 

X — , 3 

with q, = —Q^, the momentum transfer as above. In these terms, the cross section (again in 
the nucleon rest frame) is analogous to the one for elastic scattering. 
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where S is the overall center of mass energy squared. The factor in square brackets is again 
proportional to the cross section for scattering from a point charge by single-photon exchange. 
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The complexities due to the hadron target are all in the structure functions. The structure func- 
tions, however, turn out to be simpler than might have been expected. To a good approximation 
in the original experiments, the F's were found to be indpendent of momentum transfer. This 
surprising result, called scaling, was in contrast to the rapid decrease of elastic form factors. 



Doesn 7 matter much. 
Same as: 




Figure 1: Picture of deep- inelastic scattering 

It is as if, for the purposes of inclusive scattering, the Q-dependence is given entirely by the 
elastic scattering of an electron from a point charge. Representing these point charges by -k as in 
the figure, we can write schematically, 

aeN-.cX (x, g) ~ a^^"::^ (x, Q) X {x) , (5) 

where Tjq is a linear combination of the structure functions and where cri^^-!!fe* •* (x, Q) is the lowest- 
order cross section for the scattering of an electron on an elementary ^, with x given by (3). The 
precise cross section depends on the spin of the *'s, which are identified as the charged partons 
inside the proton [10]. The variable x can then be given a nice physical interpretation, if we 
think of the scattering in a frame where the nucleon and all its partons are moving in the same 
direction at velocities approaching the speed of light, the eN center of mass frame for example. 
In this frame, if the nucleon momentum is Pn, the parton's momentum should be some fraction 
of that, call it C,Pn, with < ^ < 1. Then if the scattered parton -k is to be on-shell, and and if 
we neglect its mass in the kinematics of the hard scattering, we find 

{^PN + qf = ^ ^ = 7^^ = x. (6) 

2pN ■ q 

Thus the scaling variable x observed in a specific event is the momentum fraction of the parton 
from which the electron has scattered. 
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We can now interpret the x-dependent structure functions for DIS as probabilities to find 
partons with definite momentum fractions, 



~ (7, 



e*^''^''^ {x, Q) X (probabihty to a find parton) . 



(7) 



The precise form of the parton-electron cross section depends on the spin of the parton through 
the structure functions. In this way, early experiments could already identify the spin as 1/2, 
and the partons quickly became identified with quarks. 

The parton model gave an appealing physical picture of the early deep-inclastic scattering 
experiments, but it seemed to contradict other experimental facts. If the cross section can be 
computed with lowest-order electromagnetic scattering, then partons must in some sense be 
weakly bound. But then what happened to the strong force? Already at that time searches for 
free quarks had come up empty, and the concept of quark confinement was taken as yet another 
sign of the strength of nuclcon's binding. How was it possible to describe such a force, strong 
enough to confine yet weak enough for scaling? 

2.2 Quantum chromo dynamics and asymptotic freedom 

Scaling posed the question to which QCD was the answer. The essential feature of QCD is 
asymptotic freedom [11], according to which its coupling, Qg decreases toward short distances and 
increases toward longer distances and times. Asymptotic freedom matches qualitatively with the 
requirements of approximate scaling, and its converse, that the coupling increases toward longer 
distances and times, is consistent with the behavior that gave the strong force its name. As we 
shall see, asymptotic freedom also provides systematic and quantitative predictions for corrections 
to scaling. Let's recall first how the strength of a coupling varies with distance scale. 
The (classical) Lagrangian of QCD is [12] 



in terms of n/ flavors of quark flelds and SU{3) gluon flelds A*^ — ^^T^, with the 

generators of SU{3) in the fundamental representation. The field strengths -P)ji/[^] = dfj,A^ — 
di^A^ + igs[A^, Ay\ specify the familiar three- and four-point gluon couplings of nonabelian gauge 
theory [13]. We will not need their explicit forms here. This Lagrangian, with its exact color 
gauge symmetries, also inherited the approximate flavor symmetries of the quark models, and 
was automatically consistent with properties and consequences of hadronic weak interactions 
that follow from an analysis of vector and axial vector currents [14]. In addition to the right 
currents, QCD had just the right kind of forces. 

A simple picture (if not particularly convenient for calculation) of the scale dependence of 
the coupling is shown in Fig. 2, which represents the perturbative series for the amplitude for 
fermion- vector coupling. We imagine computing this quantity in coordinate space, which amounts 
to integrating over the positions of vertices. At lowest order, this amplitude is just the coupling 
(e in QED, Qs for QCD, etc.). At the one-loop level, the amplitude requires renormalization 



/=i 



(8) 
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because whenever two vertices, say at x and y, coincide, a propagator connecting them diverges. 
For example, the coordinate space propagator for a massless boson is 



^-^H^-y) ^_ . (9) 



(27r)2 {x - yY -ie J {2ttY k^ + ie' 

Alternatively, we can just notice that the integrals over the positions of the vertices are dimen- 
sionless in gauge theories. Such integrals have no choice but to diverge for ultraviolet momenta, 
where masses are negligible. 

In these terms, the renormalized coupling may be defined (in Euclidean space) as the sum of 
all diagrams with vertices contained within a four-dimensional sphere of diameter cT, with c the 
speed of light and T any fixed time. Let's call the coupling defined this way gs{h/T), where h/T 
is the corresponding energy scale, which is then the renormalization scale. Although we cannot 
compute Qsih/T) directly, we can compute the right-hand side of the equation 

T^9s{hlT) = -^^gl{hlT) + ... , (10) 

where the term shown on the right is the finite result found when one of the vertices is fixed at 
the edge of the sphere. (We shall neglect higher orders, assuming that Qs is small.) In effect, 
the divergent part of the integrals is independent of T. The constant 60 ? of course, depends on 
the theory. For QCD, we famously have ho = 11 — 2?2j/3, with nj the number of quark flavors. 
Qsih/T) therefore decreases as T decreases, and QCD is asymptotically free. The mechanism of 



cT 




Figure 2: The coupling normalized in position space 

asymptotic freedom in gauge theories is often compared to the antiscreening, or strengthening, 
of applied magnetic fields in paramagnetic materials (whose internal magnetic moments line up 
with the external field). This behavior, associated with the self-coupling of the gluons, gives 
rise to the 11 term in h^. The — 2r2//3 term is the competing influence oi Uf flavors of quarks, 
whose effect is to screen color, in analogy to screening in diamagnetic materials (whose internal 
moments oppose an applied fleld). 
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The solution to Eq. (10), expressed in terms of the QCD analog of the fine-structure constant, 
q;s(/x) = gl{iJi)/AT^ is central to the interpretation of QCD, 



l + 6o^ln(/x'V/^2) 
47r 

(11) 
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In the first expression we relate the coupling defined by time T — h/ji to the coupling defined 
by another time T' = h/fi'. In the second expression, we use the observation that in nature it 
does not matter which starting scale /x we choose; all must give the same answer for as{^'). The 
scale Aqcd = A*exp[27r/6oQ;s(A*)] is an invariant, independent of /x, and serves to set consistent 
boundary conditions for the solution of (10). 

A simple but important observation is that although the coupling gs{^J) changes with the 
renormalization scale no physical quantity, n(Q) depends on where Q stands for any 
external momentum or physical mass scale(s). On the other hand, assuming that we can compute 
n(Q) in perturbation theory it will be an expansion in Q;s(/x), and the renormalization scale will 
also appear in ratios Q/ with physical scales. For such a perturbative expansion, we have the 
consistency condition 

dU{Q,Q/^,as{^i)) ^ ^ 
d/j, 

We will come back to this relation several times in what follows. 

The qualitative picture for DIS in Fig. 1 now has a natural explanation at large energy and 
momentum transfer. During the short time that it takes the electron to exchange a virtual 
photon, the strong force acts as though it were weak, and the electromagnetic scattering of a 
quark proceeds as though the strong force were almost irrelevant. The scattered quark, however, 
starts to move away from the target, and over a time scale of the order of the nucleon size, the 
strong force acts strongly, and produces the details of the inelastic final state, as the disturbed 
system of partons reassembles into hadrons. By this time, however, the electron is long gone, and 
the distributions we measure are detemined primarily by the original, electromagnetic scattering, 
which measures, as indicated in Eq. (5), the probability of finding quarks of various momenta 
in the parton. As the process just described proceeds, the value of relevant coupling changes 
dramatically in magnitude, from small to large. 

The discovery of asymptotic freedom opened a new chapter, not only in the strong inter- 
actions, but in relativistic quantum field theory. For the first time, effects at all orders in 
perturbation theory gave rise to observable consequences, and even the qualitative features of 
experiments could not be understood without them. 

The explanation of scaling by asymptotic freedom played the same role for hadronic physics 
that the explanation of elliptical orbits from the law of gravitation played for celestial mechan- 
ics. In both cases, most subsequent applications were (and are) to much more complicated 
experiments (scale breaking, hadron-hadron scattering for QCD, the many-body problem for 
gravitation). The striking initial success led to an essentially open-ended process of learning 
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how to use the new theory. For QCD, the problem was, and still is, to calculate in a theory 
that includes strong as well as weak couphng, and whose short distance, fundmental degrees of 
freedom (quarks and gluons) do not appear as asymptotic states (where we see only hadrons). 

2.3 Learning to calculate in the new theory 

One ultimate goal of studies of quantum chromodynamics might be to explain nuclear physics, 
say at the level that quantum electrodynamics explains chemistry. But when QCD was new it 
was natural to wonder whether we could even study the elementary quanta of the theory, the 
quarks out of which we construct the electroweak currents, and the gluons that give rise to the 
forces. As a practical matter, these particles are confined in bound states characterized by strong 
forces. Although the parton model affords a lively and compelling picture of a specific process, 
deep-inelastic scattering, it was by no means obvious that DIS was not a special case, or that 
QCD had predictions that would be realized in other experiments. 

The total cross section for deep- inelastic scattering is special because it can be related directly 
to the expectation value of a product of operators at short distances. It works this way. DIS 
is initiated by the exchange of a virtual photon emitted by an electron. The photon couples 
to quarks via the electromagnetic current, J^{x), and the amplitude to produce some hadronic 
final state \X) from nucleon state is governed by the matrix elements (X| J,,(0)|pAr). In 
these terms, the impolarized, inclusive deep-inelastic cross section at momentum transfer q is 
determined by the squared hadronic matrix elements (here simplified to a trace) 

E I (x| j,(o)|p^) 1^ 5\vx-PN - 5) = / I ny) MO) \p^) . (13) 

Scaling can be understood as a property of these matrix elements [15]. The large = —(f 
behavior of the Fourier transform to momentum space is related to the short-distance behavior 
of the expectation value of the product of currents ^. These expectation values are observables, 
and we can therefore use Eq. (12) above to derive 

(^N\J'{y)JM\PN) = {y'Y' c(^^,p^-y,as{i^')j 

- {y'y'E^n{as{W)){PN-y)\ (14) 

n 

where we have exhibited the overall dimensional behavior as a factor of and where in the 
last line we have expanded in powers of the dimensionless variable pn ■ y, Si procedure known as 
the light-cone expansion [16, 17]. After the Fourier transform, the expansion m. p^ ■ y generates 
dependence on the scaling variable, x as in the structure functions of Eq. (5). The coefficients 
Cn, however, have the interpretation of renormalized expectation values of local operators. The 

^More specifically, it is the small-y^ or "light-cone" behavior. 
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y'^, and hence g-dependence remains calculable. This is an exceptional case, however, and this 
analysis depends on our ability to sum over states in Eq. (13). 

The perturbation theory-friendly structure of the fully inclusive cross section is to be con- 
trasted to a generic S-matrix element, which depends on a variety of kinematic variables Qi and 
mass scales rrii. It will not generally be possible to pick a single renormalization scale to absorb 
the dynamical dependence of more than one such scale, and perturbation theory has little or 
nothing to say about most S-matrix elements in QCD. 

For some time after the discovery of asymptotic freedom, it was unclear whether its use would 
be limited to those relatively few, fully inclusive processes that could be reduced to expectation 
values of products of currents, or whether perturbation theory could ever say anything useful 
about the structure of final states. Correspondingly, it was unclear whether the quarks and 
gluons of QCD would manifest themselves in experiment at all. The search for footprints of 
partons took at least part of its inspiration indirectly, from an unexpected source. 

2.4 The structure of final states: cosmic rays to quark pairs 

In the mid 1950s, not long after the time when nonabelian gauge theories were first developed by 
Yang and Mills [13], jets began to be observed in cosmic ray events [18]. The term was applied 
to sprays of particles emerging from the collisions of energetic cosmic rays in an emulsion target. 
The particles typically had large (multi-GeV) and nearly parallel longitudinal momenta, with 
relatively small transverse momenta, typically of order 0.5 GeV. Such events are still difficult 
to describe fully. In the context of the parton model, however, their very structure suggests 
that, even if individual quarks might not be observable, they might still give rise to jet-like 
configurations of hadrons. 

The development of a GcV-range electron-positron collider at SLAC made it possible to test 
this idea directly. In the quark-parton model, even without QCD, a pair of quarks qq could be 
produced by the annihilation of an electron-positron pair through a virtual photon, 

e+e- ^Y{Q)^q + q- (15) 

Now the total cross section for e'^e~ annihilation into hadrons can be put into a form very much 
like Eq. (13) for deep- inelastic scattering, in terms of the products of currents. As a result, the 
total cross section is one of those special quantities for which perturbation theory can be applied 
directly. On the other hand, the lowest-order diagram for the processs (15) gives the prediction 
that the quarks will emerge from the scattering with a 1 -|- cos^ 9 angular distribution relative 
to the beam axes in the center of mass frame, with 9 the angle between the incoming electron 
and the outgoing quark. If the quark and antiquark each give rise to a jet of particles in their 
respective directions, the structure of final states will refiect directly the underlying partonic 
process. To represent the process in the parton model it was necessary to implement a model 
of hadron production within the jet that included a cutoff in transverse momentum [19]. The 
picture of the transition from partons to hadrons is represented by Fig. 3. 

Would this picture be realized? It was, once the center of mass energy reached about 7 GcV 
[20] , although to see it at these energies took some analysis. Energies have increased a hundred- 
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Figure 3: The parton to hadron transition for two jet events. 

fold in the intervening time, and jets have become a commonplace in high energy physics, in the 
final states of e"'"e~ annnihilation, of deep-inelastic scattering, and of hadron-hadron scattering, as 
illustrated by the next three figures. In deep-inelastic scattering, for example. Fig. 4, the so-called 
"current jet" associated with the scattered quark is clearly visible. Fig. 5 shows how well-defined 
two-jet events can be in electron-positron annihilation at energies around one hundred GeV, and 
Fig. 6 shows a nice example of jets from proton-antiproton scattering in the TeV energy range. 
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Figure 4: A current jet in DIS. From the HI Collaboration. 
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Figure 5: Two jet event in e"'"e annihilation. From the OPAL Collaboration. 



2.5 Jets in QCD 

By the time jets were first observed, asymptotic freedom had been discovered, and the basics of 
the analysis of DIS, which we will review below, had been developed. The parton model had 
predicted the angular distributions for jets that are initiated by electroweak processes such as 
e"'"e~ annihilation and deep-inelastic scattering, but gave little guidance on what those jets would 
actually look like, although it was thought that, as in cosmic ray events, transverse momenta 
might be limited on the scale of some fraction of a GeV. At first, it was unclear what QCD 
would have to say on this matter, until it was realized that perturbative field theory can make 
well-defined predictions for certain observable quantities even though they cannot be written in 
terms of correlation functions at short distances. These observables are not fully inclusive in the 
hadronic final state, but they do sum over many final states. They are called "infrared safe" 
[21, 22]. Suitably-defined jet cross sections are part of this class. To motivate infrared safety in 
QCD, we briefly recall the treatment of infrared divergences in quantum electrodynamics. 

In QED, infrared divergences are a direct result of the photon's zero mass. One-loop correc- 
tions to QED processes are typically finite if the photon is given a small mass A, but diverge 
logarithmically as A is taken to zero, through terms like 

Acri^l,^ ((5,me,m^ = A,aEM) ~ "EM o-i°ls((5,me)/3AB(<5Me) In^ . (16) 

In this expression, (Ja^b is a lowest order cross section, characterized by momentum transfer Q 
and (3ab is a function that is independent of A. The well-known solution of this problem [23] 
is to introduce an energy resolution, AE ^ Q, and to calculate cross sections that sum over 
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Figure 6: A "lego plot" representation of dijets in proton- antiproton scattering, represented in the 
space of rapidity and azimuthal angle with respect to the beam axis. From the DO Collaboration. 

the emission of all soft photons with energies up to AE. This is a realistic procedure, since any 
measuring apparatus has a minimum energy sensitivity, below which soft photons will be missed 
in any case. The outcome of this procedure is that corrections to cross sections including an 
energy resolution take the form, 

AE 

AaAB {Q, rUe, AE, q;em) ~ "em (^aLb{Q, ^e) f^ABiQ/me) In . (17) 

The logarithm of the photon mass has been replaced by a logarithm of the energy resolution. 
Notice that if the ratio AE/Q is not very very small, and if the function (5ab is not very large, 
ttEM ~ 1/137 ensures that (17) is numerically small. This explains the success of many lowest- 
order calculations in QED, in the face of that theory's infrared divergences. 

The Bloch-Nordseick procedure in QED suggests that something similar might be possible in 
QCD. In this case, however, since we want to use perturbation theory at high energy, we must be 
able to calculate in the limit where the ratios of light quark masses to the energy scale become 
vanishingly small, along with those involving the gluon mass, which is zero to start with. This 
leads to a whole new set of divergences, called coUinear divergences, which we will have occasion 
to discuss at length below. Something similar happens if we take the electron mass to zero in 
Eq. (17), in which case the functions /3ab develop logarithmic singularities. 

Neglecting for the moment heavy quarks, we are after a set of cross sections that are per- 
turbatively finite even when all masses in the theory are set to zero. To go to high energy in 
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perturbation theory, we must thus check that we can go to the zero-mass hmit. Once we can 
identify a quantity with a finite zero-mass hmit, and have traded zero mass back for high energy, 
we have a situation that is perfect for QCD. We will be able to use (12) to pick the coupling 
at the scale of the energy, and asymptotic freedom will ensure that as the energy scale grows, 
the relevant coupling will decrease. Perturbative predictions will then improve with increasing 
energy. 

The classic analyses of Kinoshita and of Lee and Nauenberg [24] showed that total transition 
rates remain finite in fully massless theories because the zero-mass limit does not violate unitarity 
in perturbation theory. Infrared safe cross sections are generalizations of this analysis to less 
inclusive observables. For QED, this can be done with an energy resolution; for QCD in the 
zero-mass limit, this is not sufficient. For e"'"e~ annihilation, however, we can identify infrared 
safe quantities by introducing an additional resolution. The motivation is completely analogous to 
the QED case. In the limit of zero quark mass, a quark of momentum p, = can emit an gluon 
of momentum xp, < a; < 1, [xp^ = and remain on-shell, since the remaining momentum 
(1 — x)p is still lightlike with positive energy. The resulting quark and gluon, however, are 
exactly collinear in direction, and it is by no means clear how to resolve them, especially since 
the emission, or its inverse, can take place at any time, even within a hypothetical detector. The 
same would be true for a massless electron and collinear photon. 

If we draw an analogy to the energy resolution of QED, we are naturally led to seek observables 
with angular as well as energy resolutions for high energy QCD (or massless QED), as represented 
in Fig. 7, where the cones show an angular range into which large energy fiows, while the small 
ball in the remaining directions represents an energy resolution. Without going into detail yet. 



such cross sections are infrared safe, and depend only on the overall energy Q, the angular 
resolution 6, and the energy resolution eQ, with e a small but finite number. Because they are 
physical quantities, the perturbative expansions for the corresponding cross sections satisfy Eq. 
(12), and we can write 




Figure 7: Cone jets for e"'"e annihilation. 




e,as(<5)) 
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= Y,^n{5,e)a:{Q). (18) 

n>0 

We can thus calculate these cross sections directly in perturbation theory [21], without needing 
to impose cutoffs in transverse momentum. Infrared safety alone is sufficient to produce jet 
structure at high energies, since each emission outside the cone is suppressed by a factor of 
ttsiQ). The two-jet cross sections of e"'"e~ of Fig. 7 can be generalized in a number of ways that 
we will discuss below. 

We may formalize the concept of infrared safety in the following way. QCD perturbation 
theory gives self-consistent predictions for a quantity C when C: 

• is dominated by short-distance dynamics in the infrared-regulated theory; 

• remains finite when the regulation is taken away. 

We will review in a little while how dimensional regularization can be used to control infrared, 
as well as ultraviolet, divergences. 

In "practical" terms, any perturbative cross section depends on all the mass scales (col- 
lectively, m) of the theory in addition to the relevant kinematic variables (collectively, Q) 
C{m/ fi,Q/ fi,as{fJ')). An infrared safe quantity is one whose perturbative expansion gives fi- 
nite coefficients in the limit of zero fixed masses, m, up to power corrections in m/Q, 

Cim/^i, Q/^i, «,(//)) = Yl cn(0, Q/i^) a^M + O 

n>0 

with p > 0, typically an integer. Depending on the energy range and the quantity in question, 
the mass of a heavy quark, one well above the QCD scale, may be treated as an m or as a Q. 

3 Introduction to IR Analysis to All Orders 

To explore further the theoretical basis of perturbative analysis we will work toward the proof 
of infrared safety for certain cross sections, including those for jets in e"'"e~ annihilation. Besides 
providing results that are interesting in their own right, this will also provide methods for the 
analysis of electron-nucleon and nucleon-nucleon scattering, as well as for processes involving 
heavy quarks. 

To be specific, we will develop the following necessary ingredients, sketching, wherever space 
allows, arguments that are applicable to all orders in perturbation theory. 

• A review of the treatment of UV and IR divergences in dimensional regularization; 

• A method to identify infrared sensitivity in perturbation theory: "physical pictures" ; 

• A method to identify IR finiteness/divergence: "infrared power counting"; 

• An analysis of cancellation in cross sections: "generalized unitarity" . 

The purpose of these discussions will be to explain the basis of, not to perform, explicit 
calculations. 



m 



(19) 
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3.1 UV and IR divergences in dimensional regularization 
3.1.1 Navigating the n-plane 

We will review a low order example of the use of dimensional regularization in due course, but 
even before that we should clarify how this technique can be used to treat both the ultraviolet and 
infrared sectors of QCD. The list below outlines in short the logical path from the Lagrangian of 
QCD to physically observable quantities using dimensional regularization. It also helps to explain 
why infrared safety is essential in the comparison of perturbative calculations to experiment. 



JO.QCB ^ G^''^\p,,...Pn), Re(n)<4 

^ G^'^'''\pu...Pn), Re(n)<4 + A 

^ giunphys) (p^ , . . . p^) , 4 < Re(n) < 4 + A 

^ r^'''''\{p^}) , n = A. (20) 

First, by means of path integrals or the interaction picture, the Lagrangian £qcd generates 
a set of rules for perturbation theory, out of which we may construct the Fourier transforms of 
time-ordered products of fields, the Green functions, G^'^^^^pi) as perturbative integrals over loop 
momenta. At this stage, we regularize the theory by evaluating the loop integrals in a reduced 
number of dimensions, n < 4, and by keeping all momenta away from particle poles {pf ^ m^). 
If we do this (keeping away from points where subsums of momenta vanish), the only poles in 
the n plane are ultraviolet poles. These off-shell, regularized Green functions are defined only for 
n < 4. In these terms, renormalization is a systematic procedure for removing all singularities 
at n = 4 in regularized Green functions. The resulting renormahzed Green functions, G^'^^^\pi) 
may then be analytically continued to any n whose real part is less than or equal to 4 -|- A, for 
some A > 0. 

For an infrared finite theory, this would be enough, and we could safely set n = 4, and start 
evaluating S-matrix elements and other physical quantities. For QCD, however, our job is not 
done yet, because at n = 4 the Green functions encounter a host of infrared, including coUinear, 
singularities. As a result, its perturbative S-matrix elements are hopelessly ill-defined. Our next 
step is rather to continue into the window 4 < Re(n) < 4 + A. In this range, we are still to the left 
of the ultraviolet divergences of the renormahzed theory, and (as we shall see below) to the right 
of infrared divergences, which appear in dimensional regularization as poles with Re(n) < 4. In 
this infrared-regulated range we can actually take the limits — > ml and construct the S-matrix, 
for all quarks and even for the massless gluon. 

The S-matrix elements computed for Re(n) > 4 are not physical of course, because they have 
been computed in QCD in more than four dimensions. They are not even "close" to the physical 
quantities of QCD, as is manifest by their ubiquitous l/(n — 4) poles. By contrast, however, if we 
can systematically construct quantities T(n), which are at the same time observable and free of 
poles at n = 4, the r's can be returned smoothly to the physical theory. The criterion that they 
lack infrared poles may be taken as an indication that they are dominated by the short-distance 
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behavior of the theory. These are the quantities to which asymptotic freedom can be apphed to 
provide infrared safe predictions that we may compare to experiment. Let us see how these ideas 
are reahzed for the simplest of cases, the total cross section for e"''e~ annihilation to hadrons. 



3.1.2 The total annihilation cross section 

The process e^e~ —>■ hadrons already proceeds at zeroth order in as through annihilation to 
quarks, via the photon and the Z: e"^e~ — > quark + antiquark, as in Fig. 8. Specializing for 




Figure 8: Lowest order e^e annihilation to quarks. 



the moment to a virtual photon, we can readily derive from Fig. 8 the total cross section for 
unpolarized annihilation at center of mass energy Q, 

^(0) - AT ^^'^EW n2 foT\ 

^total - ^Vcolors 2^ 

which is, in fact, a pretty good approximation to the total cross section below the Z pole, if 
we count in the sum over flavors those quark pairs that can be produced at a given energy. In 
this way, Eq. (21) is a nice test of quark masses and charges, Qjr = 4/9, 1/9, and A^'coiors = 3. 
From our previous discussions, we expect this lowest-order perturbative prediction to make sense 
for the same reason that lowest-order perturbative cross sections make sense in QED. The cross 
section cxtotai is infrared safe, and as a result, corrections will appear as a power series in as{Q), 
which is relatively small once the total energy is well into the GeV range. Let's see how this 
works at first order in a^. 

Diagrams relevant the order-a^ corrections to the zeroth-order cross section are shown in Fig. 
9. The first QCD correction (next-to-leading order, or NLO) is 

total ~ ° total ^ ° total ! l^^J 

the sum of contributions from both three-particle (quark-antiquark-gluon) and two-particle (quark- 
antiquark) final states. 

Consider first the three-particle final state. Its contribution to the cross section is proportional 
to the lowest order cross section (21), 

"i^'-^Z. (^) (^)h(Q), (23) 



TT / \ TT 
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Figure 9: Diagrams for the order Us correction to the annihilation cross section. 

where Cp = {N^ — l)/2A^c = 4/3 for QCD, and Is^Q) is proportional to the integral over three- 
particle phase space of the squared amplitude for gluon emission (here in four dimensions), 



hiQ) = (2vr) 



Q/2 



dkk d6 sin 6 
Jo 



(Q - 2k) 



X 



{Q-k{l-u))\ 
[Q-{l-u)kf ^ Qjl + u) 
A;2(l-n2) 



iQ-2k){l-u) 



(24) 



Here k is the gluon energy, 6 = Op^^ is the angle between the quark and gluon momenta, and 
u = cos 6. In this form, collinear divergences are manifest at m — ^ ±1, and infrared (soft gluon) 
divergences at A; — > 0. Notice that in this case there is no collinear divergence at k = Q, which 
corresponds to the quark and antiquark becoming parallel, or at zero momentum for the quark 
or antiquark. 

To exhibit the leading divergences, we simplify (24) to an approximation relevant to small k 
and 9: 



, , ( dk f d9 



(25) 



We now briefly review how dimensional regularization takes care of the singularities of these 
integrals. 



3.1.3 Infrared dimensional regularization in a nutshell 

We can develop the essentials of dimensional regularization as applied to infrared divergences by 
the following steps. 



'3,IR 



(2vr) 



dk r dO 
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where in the second hne we reexpressed the factor 27r as the integral over the azimuthal angle. 
To motivate the third expression we simply observe that the azimuthal angular integral is just 
a special case fli of flm, the angular integration in polar coordinates for m + 1 dimensions, and 
that the volume element in m + 2 dimensions is precisely (A; sin for m any integer ^ 

With this in mind, we reinterpret the parameter m in Eq. (26), which in the unregularized 
theory has the value unity, as a complex parameter, and redefine flm as a complex function fl{m) 
of m, restricted to be the angular volume ftm in m + 1 dimensions for any integer value of m. 
This turns out to be enough to define Jl(m) uniquely for all complex m. 



27r(rf+i)/2 



r(i(rf + : 

= 2V/^^' (27) 

with r{z) the Euler Gamma function, with the well-known properties, 

r(i/2) = v^, 
^r(^) = T{z + i), 

r{N) = (iV-l)!, (28) 

where N is an integer. With these properties we easily check the angular volumes for low 
dimensions: m + 1 = 2 {27i), 3 {An), 4 {27i'^). 

For Re(m) > 1, /3.1R becomes finite, and the soft and coUinear divergences appear as poles 
at m = 1, as anticipated. 

This proceedure defines a new theory: QCD in 4 + (m — 1) dimensions, where m = 1 is real 
QCD. 

A similar continuation in dimension can always be found for any observable, since whenever 
we continue expressions to n > 4, we will be able to integrate over the angular phase space 
for dimensions above n — A. This, in turn, will regulate the remaining integrals. Analogous 
considerations apply to the ultraviolet divergences regulated with n < 4 as in Eq. (20) above 
[25]^ 

^This is easy to verify inductively, using polar coordinates in N dimensions to define cylindrical coordinates 
in A'^ + 1 dimensions as an intermediate step. 
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Our notation below will generally be to take n as the number of dimensions, and to express 
variations away from n — A through e — 2 — n/2. Again, we emphasize that it is only infrared 
safe quantities that can be evaluated at £ = 0, where they return to the "real" theory. 

For the case at hand, in n = 4 — 2£-dimensional QCD, wc found double and single poles in 
/a JR. The full expression is also straightforward to evaluate, up to terms that vanish for £ — > 0, 
and cr^otaf^'* given in this approximation by 



Clearly, the three-particle contribution to the total cross section is highly singular when returned 
to four dimensions. Although unphysical, it is positive, corresponding to the absolute square of 
the amplitude for the emission of a single gluon. 

Dimensional regularization, however, also enables us to compute the two-particle contribution, 
given by one-loop vertex and self-energy diagrams. This calculation is perhaps more familiar, 
since it can be carried out in terms n-dimensional loop integrals. Unlike the three-particle cross 
section, it is not an absolute value squared, and need not be positive. The result is proportional 
to the lowest order cross section, and is given by 



Here wc sec the same poles as in the three-particle case, but now with negative rather than 
positive signs. As expected the sum of the two- and three-particle NLO cross sections is finite, 
and the total NLO cross section is very simple. 



a satisfying result that is nonsingular £, and can hence be returned smoothly to four dimensions. 

Our goal is now to identify classes of cross sections, involving jets and related observables, 
for which we know ahead of time that calculations like these will provide predictions at £ = 0, 
that is, for real QCD. To do so, wc need first to develop a better understanding of where poles 
in £ come from, not only at NLO, but to all orders in the coupling. 

3.2 Physical pictures 

Divergences in perturbation theory can come about in only two ways: cither from the infinite 
volume of momentum space, or from infinities in the integrands of individual diagrams. In 
renormalized perturbation theory, the former are eliminated by the systematic construction of 
counterterms. The latter, of course, reflect the vanishing of the causal propagators: k'^ —m^+ie — 
0, where m may or may not be zero. We have kept the infinitesimal imaginary part of the 
propagator to remind ourselves that perturbation theory is defined by integrals in the complex 
plane of every loop momentum component. The integrand of any diagram is thus a meromorphic 




(30) 




(31) 




(32) 
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function of the 4L complex variables of integration in four dimensions, where L is the number of 
independent loop momenta. This will be a very useful observation, as we shall soon see. 

To begin with a specific example, consider the one-loop quark electromagnetic form factor, 
the sum of the lowest-order vertex and self-energy corrections to the electromagnetic current, 

^) = -e/^^ u {pih„v{p2) e) . (33) 

The function e) is the correction that appears times the lowest order annihilation cross 
section, Eq. (31) at NLO for the two-particle final state, 

M^'^) - -^^^[-Z^e) r(l - 2e) IF^" 2(^ + 7- 

Because the electromagnetic current is conserved, the ultraviolet divergences in the sum of dia- 
grams cancel, and no counterterms are necessary in this computation. 

The infrared poles of dimensional regularization, which diverge at n = 4, require, as pointed 
out above, that some lines go on-shell lines for some real values of the loop momenta. But on- 
shell lines at real momenta are not enough. The momentum integrations, interpreted as contours 
in complex planes must also be pinched between coalescing poles, at least one pole on either 
side of the contour. If not, Cauchy's theorem can be used to deform the contour(s) away from 
the poles so that the integral can carried out in such a way that the integrand is always finite. 
The requirement that the contours be pinched makes it possible to identify potential sources of 
infrared singularities systematically. 

Consider a particular subspace of the momentum space of a given diagram where some set of 
propagators diverge, that is, where some set of lines are on-shell. We'll refer to such a subspace 
as a "singular surface" of loop momentum space. Let us contract all the fines that are not on- 
shell at the singular surface to a point. The resulting diagram is called a reduced diagram. The 
singular surface can give rise to a divergence only if all the loop momenta of the reduced diagram 
are pinched between coalescing singularities at the surface. We therefore refer to such a surface 
as a "pinch surface" . The existence of a pinch surface is a necessary condition for the production 
of an infrared pole. 

It is natural to think of the reduced diagram, all of whose lines are on shell as the picture of a 
process, but in fact this is not automatic. In a physical process, particles travel between points, 
represented by the vertices where they begin and end. In these terms, we can use a criterion 
identified by Coleman and Norton [26, 27] and distinguish pinch surfaces solely on the basis of 
their reduced diagrams: the lines of the reduced diagram must describe a scattering process in 
which classical point particles move between vertices, such that all their motions are mutually 
consistent. 

To be specific, for a reduced diagram with lines of momenta pk to correspond to a pinch 
surface, we must be able to assign to each of its vertices i a position Xi in space-time such that 
if vertices i and j are connected by line k, with x° > then 

(x,-x,r = (x°-x;)<, (35) 
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where Vk = Pk/Pk the four-velocity associated with the on-sheU momentum p^. This corresponds 
to free propagation between vertices in space-time. The condition (35) is invariant under a scale 
transformations for all points, keeping momenta fixed: Xi — > Xxi, for all i. For arbitrary A > 1, 
the process simply "gets larger" , and we can think of this as a potential source of the infrared 
divergences. 

The proof of this assertion is easy, and may be demonstrated adequately by a simple example, 
the triangle diagram in Fig. 9. Because the divergences are associated with momentum depen- 
dence in denominators, we simplify further by suppressing all numerator factors, and just study 
the scalar triangle diagram with light-like external momenta, pi and p2, 

/a(pi,P2) = / + ie) ((pi - + ie) {{p2 + kf + ie) ' ^^^^ 

which, introducing Feynman parameters, can be expressed in terms of a single denominator, 
that is quadratic in the momenta, 

^ . . ^ [ dr-k daida2daz5{l — YH=i'^i) 

= 2y^y^ — ^ 

D = aik^ + a2ipi - kf + a^{p2 + kf + ie . (37) 

Because D is quadratic, its zeros, and therefore the points at which the integrand diverges, 
are simply the solutions to a quadratic equation in each momentum component. A necessary 
condition that each momentum component be pinched between singularities is that the two 
solutions for each k^ must coincide, = ... 3, so that 

^D{ai,k^,Pa)^0. (38) 
For the case at hand this leads to the following conditions, 

aik — a2{pi -k) + az{p2 + /c) = 

t 

AtiVk - At2Vp^-k + AtsVp^+k = , (39) 

in terms of times given by Ati = Q:ip^ and velocities Pi/p'^, as noted above. The second condition 
is the statement that if wc start at any one of the three vertices, the total four-distance travelled 
by going around the loop is zero. This is equivalent to the statement that the points are connected 
by classical propagation for each of the three particles. Now in fact, for specific solutions, one or 
two of the times may vanish. Indeed, whenever a line is off-shell, its corresponding a parameter 
must vanish, or we may deform the a contour and avoid the point D — without worrying at 
all about the momenta. The points ai = 0, however, are boundaries of the a integrals. These 
boundaries define the integrals, and therefore cannot be avoided by contour deformation. The 
vanishing of an ai normally corresponds to the vanishing of the time of propagation for an off-shell 
particle, one contracted to a point in the corresponding reduced diagram. 

There are three solutions for Eq. (38), which illustrate the generic physical interpretations of 
e+e" annihilation processes. 
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• CoUinear to pi : in which two particles of momentum pi — k and k are on-shell and exactly 
parallel to one of the external lines. The third line, with momentum (p2 + k) is off-shell, 
and its time of propagation is correspondingly zero 

k = CPi, ct3 = 0, ctiC = a2{l - C) ■ (40) 

• CoUinear to 

k = -Cp2, a2 = 0, aiC = o;3(l - CO ■ (41) 

• Soft: for which the solution is ai = 1, while the velocity of the zero-momentum particle 
is undefined. There are similar pinch surfaces when the remaining two momenta vanish. 
The physical picture in this case may be thought of as an infinite-wavelength soft particle 
coupling to finite momentum particles at arbitrary points in space-time, 

A;^ = , (aa/ai) = (ag/ai) = . (42) 

These configurations show that the sources of infrared sensitivity that we identified above for 
the high-energy and/or zero mass limits, collinear and soft, correspond directly to pinch surfaces, 
at least for the example that we have examined. The same analysis, however, can be carried out 
in just the same way for an arbitrary diagram with arbitrary numbers of loops. This leads to the 
generalization of Eq. (38) known as the Landau equations [28], 

foralH : £^ = m^, and/or ctj = 0, 

and OiiliSis = 0, (43) 

i in loop s 

where =1 (-1) for loop momentum i flowing along (opposite to) momentum £j, and otherwise. 
Generalizing Eq. (39), the consistency of the physical picture associated with an arbitrary pinch 
surface is ensured by the identity 

Axi2 + AX23 + . . . + Axni = (44) 

around each loop of the diagram, of the sort in Fig. 10. 

The interpretation of pinch surfaces, and hence of the possible sources of infrared poles, in 
terms of physical pictures greatly simplifies the process of classifying their origin in many cases, 
particularly for processes involving hard scatterings. As we have suggested above, however, 
satisfying the Landau equations is still only a necessary, and not sufficient condition for an 
infrared divergence. Clearly, the strength of the singularity on the pinch surface plays a role 
as well. This leads us to power counting, which provides us with a method for distinguishing 
singularities that can give rise to divergences from those that cannot. 
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Figure 10: Loop illustrating the physical condition (44). 
3.3 Power counting 

As with our analysis of physical pictures, our goal in analyzing the strength of singularities by 
power counting is to find necessary conditions for divergences. Until we do a specific integral 
in detail, it is always possible that an unexpected symmetry or cancellation between diagrams 
might eliminate an apparent singularity. Indeed, in gauge theories conspiracies between diagrams 
are more common that not, because only in gauge- invariant sums do unphysical modes decouple 
from physical quantities. Without going into detail right here, it is worth pointing out that 
unphysical modes often give rise to unphysical infrared poles in individual diagrams [29]. 



3.3.1 Example: the triangle with scalar quarks 

We once again turn to the triangle diagram of Fig. 9 to illustrate the method. This time, however, 
we will analyze the diagram with massless scalar "quarks" coupled to a massless "gluon" or 
photon, taken for simplicity in Feynman gauge. 



^IR 



iR [-2{pt -k+)k- -kl^ie) (2{p-^ ^k-)k-^ -kl^ie) {2k+k- -kl^ie) 



Here we are using light-cone momentum components, 

k- = ^^{k^^e) 

k± = {ki,k2), (46) 

with external momenta Pi — (5^+ , P2 — 5^- ■ 

We anticipate on the basis of pinch surface analysis that Air is singular when the gluon's 
momentum k is soft, that is, vanishing in all four components, and also when it is collinear to 
either pi or to p2. It is also clear that these regions overlap. We would like to develop a technique 
that enables us to verify the behavior of the diagram in these limits, without double counting. 
Correspondingly, we would like a means of putting bounds on the integrals in the presence of the 
singular points. We begin by having a look at the limit of soft momentum. 
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Wc can zero in on soft momenta by the following rescaling, where we'll be interested in what 
happens when A vanishes, 

n^'^Xki'. (47) 
We implement the rescaling by inserting unity in Eq. (45), in the form 



isoft^ [^^^ d\H [>? -Y^K^ . 



(48) 



The argument of the delta function is by no means unique; at present it is preferable to choose it as 
a polynomial in momentum components and therefore an analytic function of of the components. 
Inserting (48) into (45), we then find the (noncovariant) expression 

Asoft ^ _ / dX^xn / fin^ (4^+^- + o(A)) 5 (l - - - 4) 

1 

^ A(-2p^K- + 0(A) + ie) A(2p^fi;+ + 0{X) + ie) X^{2k+k- - + ie) 

. _2 / , ^ ^'w" . (49) 

io A^-" J {-K+ + le) {k- + le) {2k+k- -K^ + ie) ^ ^ 

We have dropped terms that are higher order in A in both numerator and denominator. The 
integral is then homogeneous in A, and the A integral can be performed to give a 1/s pole in n 
dimensions. The remaining integral is free of soft divergences since is fixed at unity. The soft 
divergence, where all components of momentum vanish together, has therefore been isolated by 
the scaling (47). This does not mean that the remaining integral is finite in four dimensions, 
however, and we easily verify additional singularities associated with the presence of coUinear 
divergences at finite k^, and arbitrarily close to the point — 0. These can be treated in much 
the same fashion as for the original diagram, by using the delta function to do one of the k 
integrals, and then analyzing the analytic structure of the remaining integral. 

The additional coUinear divergences appear in the k integrals as pinches at the surfaces defined 

by 

= 1, 

= 4 = 0, (50) 

where we note that the variables k are dimensionless. These are the soft tails of the coUinear 
regions, and following a similar analysis we find that they produce double poles. For example, 
for the pinch surface where = 1 we may introduce a second scaling, 

K-' = A' K~, k't? = A' 4, A'2 = + (4)2 , (51) 

and verify that the A' integral, which now controls the approach to the pinch surface, produces 
a seond pole in e. The remaining integrals are then manifestly finite. We discover in this way 
the double poles for the scalar quark vertex; going through identical steps in QCD will give the 
same result. We now turn to a more general analysis. 



24 



3.3.2 General pinch surface analysis 

As in the discussion of the Landau equations, the steps from the specific example to a general 
analysis are clear. Suppose we consider an arbitrary graph G{Q) with specified external momenta, 
denoted collectively by Q. Using the physical picture analysis, we can in principle identify 
every pinch surface 7 of G{Q). At any point on 7, we may divide the coordinates of the loop 
momentum space into "internal" coordinates, i^, that specify directions within the surface 7, and 
"normal" coordinates, that specify directions that move us off the surface. For example, for 
the soft divergence of the triangle diagram the pinch surface is a point, so that for this surface 
all corrdinates are normal coordinates. For the coUinear-pi pinch surface, with pi in the 
plus direction, the plus component of the k loop momentum is internal, while the minus and 
transverse are normal. 

Isolating the integral of G{Q) near 7 the same way we treated the triangle diagram near the 
soft pinch surface (point), we write in dimensional regularization 

where the numerator factors and denominator factors dj are polynomials in the normal and 
intrinsic coordinates, along with the external momenta. Any Jacobean factors are absorbed into 
the definitions of the volume elements, and the parameter ^^^{e) represents the number (possibly 
noninteger) of normal variables. Once again, we insert unity in the form of a scaling for the 
normal coordinates (only), Ka — A-yK^, 



i 




(53) 



As above, the choice of normal coordinates will depend on the pinch surface. 

Under the scaling, both numerators and denominators have dominant terms in the A — > 
limit, the terms of lowest overall powers and Dj of normal coordinates, respectively, 

n.(/€a,4,g) = Af [n,«,4,g) + 0(A)] 

dj{^^a,h,Q) = X^^ [dj{Kle,,Q) + 0{X)] . (54) 

In terms of the exponents Nj, and Dj, we find in this approximation an overall scaling behavior 
as we approach the singular surface, given by 

rK""" d\ 

P7 = J2Dj + l-V^{e)-J2Ni. (55) 
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A^{Q) is the remaining integral, in which every numerator and denominator factor has been 
replaced by its term of lowest homogeneity in the k's, 




J 2 



a 



) 



(56) 



On this basis, we estimate the contribution from the pinch surface 7 in dimensional regularization 
as the result of the A integral, times whatever the remaining integral A^((5) gives. 

In general, the integral A^((5) will have additional pinch surfaces, of course. If, however, these 
pinch surfaces are among the original set of surfaces identified for the graph G{Q), we can bound 
the behavior of A^{Q) on this basis, and eventually bound the entire integral by an expression 
that has some maximal number of poles (possibly none) in e — 2 — n/2. This situation is by 
no means assured, however, and different classes of diagrams require different choices of normal 
coordinates. It is by this reasoning that we verify, however, the infrared safety of classes of jet 
and related cross sections. These cross sections are in general not free of pinch surfaces, but are 
defined in such a way that these pinch surfaces are not strong enough to produce poles in s. 

We note that this analysis is simpler for massless quarks. Organizing poles in dimensional 
regularization is more complicted in the presence of quark masses, simply because the masses 
provide additional scales. 

3.3.3 Outline: general analysis for e+e" annihilation 

The ingredients above are sufficient to identify the sources of long-distance behavior in e+e^ 
annihilation with massless quarks quite generally. Pinch surface analysis begins with the obser- 
vation that all momentum fiows into the hadronic sector through a local operator, an electroweak 
current. We must identify all physical pictures where particles emerge from that point, interact 
and eventually form a (perturbative) final state. 

The reasoning is very simple, keeping in mind that all particles are taken to be massless, and 
therefore propagate at the speed of light. As a result, on-shell, finite-energy particles that emerge 
moving in different directions from the short-distance ("hard") vertex at which the current acts 
instantly become space-like separated and can never interact again to exchange finite amounts of 
momentum. On the other hand particles that emerge from the current vertex moving in exactly 
the same direction may interact, recombine or split at any time. The final state formed in this 
way will always consist of sets of exactly collinear particles for all times. This is the origin of jets. 
In addition, all finite-energy particles moving in whatever directions can always be connected by 
zero-momentum, infinite wavelength lines. Thus, although different jets of particles may not 
exchange finite momenta at a pinch surface, they can exchange color quantum numbers in a 
nonabelian theory. 

The reduced diagram for any pinch surface in e+e" annihilation then reduces to: 

• A set of off-shell-hnes at the point where the current acts. We refer to this as the hard 
scattering subdiagram (or function), H{Q). 
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• A set of jet sub diagrams, Jj of exactly collinear lines that emerge from the hard scattering 
through one or more finite-energy "jet lines" , which rearrange themselves into any number 
of particles, all moving with exactly parallel momenta. 

• A soft subdiagram, S, made up of lines whose momenta vanish at the pinch surface. The 
soft subdiagram may be connected to any of the jet subdiagrams, or to the hard subdiagram 
through zero-momentum lines. 

Power counting, carried out as above, may be extended to pinch surfaces of semi-inclusive 
cross sections, by scaling final-state mass shell delta functions in the same way as propagators, 
which have the same dimensions and homogeneity properties, <->■ l/Zc^. In this manner, 

we can show that in a gauge theory like QCD, cross sections with fixed numbers of particles, but 
with the same dimension as the total cross section are at worst logarithmically divergent, and 
have up to two poles in e per loop [27]. We will refer to such cross sections as semi-inclusive. 
They are simply cross sections for which every particle in the final state is integrated over at 
least part of its phase space. 

Space does not allow a detailed proof of logarithmic power counting, but we can think about 
it in the following way for the the quark- ant iquark production process. We have already seen 
that the one-loop vertex diagram is logarithmically divergent for n = 4. Now let's consider what 
happens when we attach an extra gluon line to the vertex diagram, so as to preserve the pinch 
surface structure. 

Adding one additional vector line to such a diagram increases the number of lines by at 
most three, and increases the number of loops by one. If the line is soft wc may choose the 
normal variables as all four of the loop momentum components. If we attach the soft line at 
each end to jet lines, the denominators of the new jet lines are linear in the soft momentum, 
while the denominator of the new soft line is quadratic. Scaling the new soft loop along with the 
original diagram we find an extra power from the new denominators and from the new 
loop momentum, so that the power counting remains unchanged. Adding a jet line leads to the 
same result, once the scaling dependence of momentum factors in the numerator are taken into 
account. 

For a process which like e'^e" annihilation has a single hard scattering in its reduced diagram, 
we may choose normal variables according to the example shown above. The reduced diagram 
for any such pinch surface will be characterized by a set of jets, which we may label by index i. 
Fig. 11 shows a typical example of a reduced diagram, with two jets for simplicity. The figure is 
presented as a cut diagram, with the perturbative contribution to the amplitude on the left, and 
the complex conjugate amplitude on the right. The pinch surface analysis and power counting 
described above is readily extended to cut diagrams through the power-counting correspondence 
of mass-shell delta functions and propagators just noted. 

The overall power of A'**^'''''^^ for pinch surface 7 may be thought of as defining a "superficial" 
degree of infrared divergence, iJ,{j) (the extra -1 is compensated by the volume element dX) . The 
precise graphical and momentum flow configuration of the reduced diagram are directly reflected 
in the value of //(7). At an arbitrary pinch surface, we may identify a number of relevant integer 
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Figure 11: Representative two-jet reduced diagrams for jet cross section in cut diagram notation: 
(a) physical gauge; (b) covariant gauge. From Ref. [30] 

parameters, associated with the hard, jet and soft subdiagrams, that provide a lower bound on 
the degree of infrared divergence. 

• Some numbers iJj > 1 of lines attach the jets in the amplitude and its complex conjugate 
(to the left and right of the final-state "cut" in a cut diagram). In the Hs we count only 
quark lines and those gluon lines with transverse polarizations at the pinch surface. (In 
covariant gauges, ghost lines should also be counted.) The contributions from pinch surfaces 
where one or more jets connect to a hard scattering by only unphysical gluon polarizations 
vanish [29] by the Ward identities of QCD, which ensure that physical states evolve into 
physical states only. 

• A\^'^^ unphysically polarized gluons may, however, appear as additional gluon lines at- 
tached to the hard scatterings. 

• The soft function S may attach to the jet subdiagram Jj, i = 1 ... J in an arbitrary 
fashion, by the exchange of soft gluons, and in the general case soft quarks and even ghosts 
in covariant gauges. We denote the number of soft gluons, quarks (antiquarks) and ghosts 
that connect S with Jj as a^, Qi and q, respectively. 
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• Some numbers of soft lines S^^'^^ may attach S to the left and right hard scattering vertices. 
(Such lines are not shown in figure.) 

The minimal superficial degree of infrared divergence for an arbitrary J-jet pinch surface in 
covariant gauges may be computed in terms of these parameters and is given in four dimensions 
by [27] 

/.(7) = + + I E {Hi'^ + h\^^ - 2 + 2g. + c. + + aS^^ + a. - «f ] } . (57) 

jets 1=1 

where u\ is the number of three-point vertices in jet i to which soft gluons Oj or jet gluons Ai 
connect. In physical gauges, the terms involving the AiS are absent. 

The minimal value of /x(7) of a given cut diagram defines the most divergent, or "leading" 
pinch surfaces (referred to as "leading regions" in Refs. [27, 30]). For negative /x(7) we have an 
overall power divergence; for ^^(7) = the overall divergence is logarithmic, and for positive //(7) 
the pinch surface does not produce a divergence at all. 

For an arbitrary pinch surface relevant to a semi-inclusive cross section in e"'"e~ annihilation, 
we easily verify that the /i(7) > 0. The lower limit, corresponding to logarithmic divergence, 
is realized only for reduced diagrams like the ones of Fig. 11. These are characterized by quite 
specific conditions on the variables above. 

• At most a single quark or physically polarized gluon connects each hard scattering vertex 
to each jet: h\^'^^ — 1. 

• Arbitrary numbers of scalar-polarized gluons, with their polarization vectors at the hard 
scattering proportional to their momenta, may attach the hard vertices to the jets in 
covariant gauges. These are absent in physical gauges. 

• Soft functions are attached to jets only by soft gluons and at three-point vertices only. 

• No soft gluons are attached directly to the hard vertex. 

Our ability to classify pinch surfaces, and to estimate their degrees of divergence, in this 
manner is the basis of effective field theories that systematically separate the dynamics of hard 
scattering, jets and soft partons [31]. 

It is worth noting that fermion and/or ghost loops are permitted internally to the soft function. 
In fact, in massless quantum electrodynamics, it is through fermion loops that photons interact 
with each other, both in the soft subdiagram, and in jet subdiagrams. 

Infrared safety requires that contributions from these well-characterized logarithmically diver- 
gent regions cancel. In the next section, we will show how this comes about in e"'"e~ annihilation. 
Our arguments, some of which are new to these lectures but which are closely related to those 
in [27, 30], will start with an analysis of energy fiow. 
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3.4 Cancellation 



As we have just seen, pinch surfaces for semi-inclusive cross sections in e^e~ annihilation are 
characterized by fixed numbers of jets. All interactions within the jet subdiagrams involve the 
splitting and/or recombining of coUinear lines, which preserve the total jet four- momenta. We 
will now show that if we sum over all final states with the same flow of energy, that is, the same 
numbers of jets, these logarithmic divergences will cancel, as long as the the sum over final states 
is not weighted by a function that is itself singular [32]. 

We begin by introducing an operator that measures energy flow in an arbitrary state with 
specified particles and momenta [33], 

n 

E{h)\ki...kn) = Y1 u;{kj)S^{kj-n)\h...kn). (58) 

To construct the operator £, we work in the interaction picture for the fields, where creation and 
annihilation operators are time-independent. Within this context, the reasoning is quite general 
for bosons and fermions, and we identify the operators by a generic index, i. In these terms, the 
energy flow operator is simply 

S{n) = Y, I dHa\{i)ai{i) S\i - n) . (59) 

We have chosen to normalize our creation and annihilation operators by 

[ai{k), a]{q)]± = 6ij2ujk6'{k - q) , (60) 

where the ± subscript indicates commutation for bosons, anticommutation for fermions. 

Using S{n), we next construct what we will refer to as an event shape operator, specified by 
an arbitrary function of angle, /(rh,n), 

:F{m) = Q-^ J d% f{m,h) £(n) = Q-^ J £1 a\(i)ai(i) f{rh,i) , (61) 

where m represents one or more fixed directions that we may decide to incorporate into the 
definition of the weight function /, and where Q has dimensions of mass. For e"'"e~ annihilation, 
we generally take Q to be the center of mass energy. 

Any state with definite particle momenta is an eigenstate of JF, with an eigenvalue given by 
the sum of the energies of each particle times the corresponding weight function, 

n ^ ^ 

T{m) \ki . . . kn) = Q''^ ^ uj{kj) f{m,kj)\ki. . .kn) 

i=i 

= fN\ki...kn). (62) 
The eigenvalue /jv is the total weight of state \N) = \ki. . .kn)- 
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Event shape operators of this kind can be used to define difi^erential cross sections in e+e 
that sum over all states of a given weight, which we may represent as 

^ = L'-(fe,fe)%^ (63) 

for the annihilation of a Icpton pair of momenta ki^2 into hadrons. The cross section is written 
as the product of a leptonic tensor, calculable in QED, and a hadronic tensor, where the indices 
are associated with the electromagnetic (more generally electroweak) currents that couple to 
the quarks. Using our operator we express the hadronic tensor as 



= E / d'xe^'^- m,{x)\N) {N\S{f - :^)>(0)|0) . (64) 

These are semi- inclusive cross sections in the sense defined above, and in general the contributions 
of states of definite particle number to these cross sections are logarithmically divergent in the 
zero-mass limit. 

Of course, Eq. (64) is rather abstract, but we can represent it more explicitly in terms of the 
energy flow operator by using the integral representation for the delta function, 

^(J - /•) = f°° ^Q-iT[f-Q-' J d^n f{m,h) £(n)] _ ^gg^ 

To generate a perturbative expansion, we use the interaction picture to evaluate the matrix 
elements in Eq. (64), 

Y: m^-'I''y^'Ux))\N) {N\5{f - ^)T(e^/'^Vr.^-^(o)) |o) 

N 

J-oo ZTT 

^ir J d^h f(m,h) £{h) T{e'I'^'y'^'j^{0))\0) , (66) 

with Cj the interaction terms of the Lagrange density. These matrix elements are evaluated in the 
free theory, according to Wick's theorem, with T denoting time ordering within the amplitude, 
and T anti-time ordering in its complex conjugate. 

At each order of perturbation theory, infrared divergences arise from vanishing phases at 
large values of the time variables i/o, y'^ for the interaction Lagrangians. This claim is adequately 
illustrated by a scalar 0^ coupling, which in terms of creation and annihilation operators gives 

d^yCi{yo,y) j d^kid^k2C{ki,k2)a\ki)a\k2)a{ki + k2) 



y^Q-iyo['^{k{)+uj{k2)-uj{ki+k2)] _j_ 

C{kuk2) = ^ 1 ^ , (67) 

8(27r)3a;(A;i)a;(A;2)cu(A;i + £cA;2) 
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where we show a representative term"^. The time yo has precisely the interpretation of the 
time at which an interaction occurs. Infrared divergences are produced only when interac- 
tions from arbitrarily large yo contribute (the integrand is just a phase, and is thus finite). 
Oscillations in the exponentials of Eq. (67) damp the integrals over yo, unless the sum over 
energies vanishes in the exponent: J2 — * 0. It is immediately clear that the combination 
—iyo[u!{ki) + uj{k2) — uj{ki + k2)] vanishes for soft and coUinear emissions only. More generally, 
the vanishing must occur at a saddle point in the integrals over momenta, corresponding to pinch 
surfaces in the corresponding diagrams [34]. 

Now, in the absence of the operator involving 8 in the middle of Eq. (66), the large yo^y'o 
dependences of the interaction Lagranians cancel. That is, if we make the replacement, 

gir / d?h f{rh,h) £{h) _^ ^ ^ ^gg-^ 

the exponents cancel for large yo,y'o- This is just the statement that the total cross section is 
free of IR poles, as required by unitarity [24]. 

There are various ways of showing this result. One instructive approach uses the optical 
theorem, which relates the total cross section to the imaginary part of the same matrix elements 
above with replacement (68), 



^%-{Q) = ^Ini7r(g2) 

7r(g^) = -(l/3)i / d'x e-^- (0| T {f{x) j^(0)) |0) . (69) 

We now invoke the pinch surface analysis to the function 7r((5^), which is the just the sum of 
photon self-energy diagrams. We recognize that there are no physical processes for these diagrams 
at any order of perturbation theory. Such a process would require a set of jets to be created at a 
point and then, without exchanging energy, to propagate finite distances in different directions 
and finally to annihilate. Such a sequence of events is simply not possible within a physical 
picture. 

We cannot do quite the same thing for weighted cross sections, just because the interaction 
Lagrangian does not commute with the weight operators in general. Nevertheless, as we have 
constructed them, the weight operators T actually do commute with the interaction Lagrangian 
at the special momentum configurations necessary for cancellation of infrared divergences. This 
works because all commute with the soft and ccoUinear parts of Ci. For example, with the 
trilinear coupling above, we can compute 

j d^kid^k2 C{ki, k2) a^iki) a^{k2)a{ki + k2) , J d^I a\{t)ai{t) f{m, 1) 
= j d^kid^k2a\ki) a\k2)a{ki + k2) 

2-'(A i + /.•2)./(m, fik.+k^) - 2u;{ki)f{rh, h) - 2uj{k2)f{m, h) 1 ■ (70) 
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^With our normalization 60 for creation and annihilation operators, our scalar field is 4>{x) = 
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The term in brackets on the right vanishes for any coUinear set of momenta, or when one of the 
momenta vanishes. This cancellation of interactions was referred to as generalized unitarity at 
the start of this section. In the same way, all event shapes functions, or weights, that are linear in 
energy times a smooth angular dependence give matrix elements in which largc-timc interactions 
cancel. Such weight functions are infrared safe. In addition, once a differential event shape is 
shown to be infrared safe, we can take averages or moments of it, extremize it with respect to 
the variable directions, m, and so on. We now turn to some examples. 

3.5 Jet cross sections and event shapes 

The most direct definition of a jet is in terms of observed energy within a cone, a choice that we 
motivated in Sec. 2 on the basis of angular resolution. Such a cross section can be implemented 
in terms of energy flow operators as above, and is infrared safe. Actually, some consideration 
has to be given to the discontinuous nature of the corresponding angular function f{rh, n) at the 
edge of the jet cone, with rh representing the cone directions and widths. This does not produce 
singularities, however, because divergences are at worst logarithmic, while the discontinuity is 
encountered as an endpoint [27]. 

Any number of variations on the theme of jet cross sections can be given. Cluster algorithms 
are a popular alternative to cone jets [35]. In a cluster algorithm, particle momenta within any 
final state are successively combined by searching for the smallest value of some variable that 
vanishes when two particles are coUinear. For example, in the "/ct" or "Durham" algorithm 
applied to e+e" annihilation^, we define i/ij for momenta Pi and pj by 



We start with the pair of particles that has the minimum value for i/ij. These two momenta 
are combined into a single momentum according to any continuous rule that amounts to simple 
addition when the angle 6ij vanishes. The algorithm is then repeated, until all remaining pairs 
of momenta satisfy the condition 



for an arbitrary, and adjustable, value of the parameter ycut- The sets of particles whose momenta 
are grouped in this way constitute the jets for that particular value of ycut- The smaller j/cut the 
more coUimated are the jets, and of course, the number of jets generally increases as ycut decreases. 
This illustrates that neither the grouping of particles into jets nor even the number of jets in a 
given final state is unique. 

Uniqueness is not necessary, however, for the concept of the number of jets to be useful. What 
is necessary is that it be possible to compute the distributions of jet number self-consistently in 
perturbation theory. An infrared safe A^-jet cross section appears at lowest order as an A^- 
parton production process. This provides a basic physical interpretation. Corrections to this 
approximation associated with different, more complex quantum histories are inevitable, but 
calculable. 

^See [36] for a fcy algorithm formulated for hadron-hadron scattering. 



Vij = 2Q-2min [{p',f, {p'.f] (1 - cos%) . 



(71) 



Vij > Vent 



(72) 
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Labeling a given state as two-jet or three-jet says almost as much about the definition of 
a jet as about the final state. It is therefore also useful to analyze final states more directly 
through weights called event shapes [37] that can take on continuous values. The thrust [38] is 

the prototype event shape. 

The thrust quantifies the deviation of a final state from a perfect configuration of two "back- 
to-back" jets, by projecting momenta on an arbitrary axis rn, and then maximizing with respect 
to m: ^ 

T = — maxrfj ^ |m ■ pi| . (73) 

When t = 1 — T approaches zero, the final state consists of two very narrow jets. Another 
interesting relation, quite easy to verify using light cone variables, is that for narrow jets, 1 — T 
is proportional to the sum of the squares of the invariant masses of the two jets, computed from 
particle momenta within the hemispheres defined by the plane normal to the thrust axis, 

t^l-T -h^r^ ■ (74) 

Equivalently, if wc define the z direction in the direction of jet Ji (chosen arbitrarily), and we 
define a standard rapidity for each particle momentum relative to that axis. 




(75) 

then, neglecting particle masses, the deviation of the thrust from unity may also be written as 



Q 



Wk,Te-M. (76) 

^ all i 



Starting from this form, we can generalize the event shape t = 1 — T to 

^a-^H kiT e-l'^l^^-"^ , (77) 

^ all i 

where eg is t — 1 — T, and ei is called the "jet broadening" [39]. As a class, the ea's have been 
named "angularities" [40, 41]. 

Figure 12 shows a variety of event shapes as measured at the LEP collider [42], including 
the thrust, and variations of the broadening. We notice that except for the two-jet limit near 
1 — T = 0, for example, the event shapes are relatively slowly varying functions. The arrows in 
each panel show the range over which the OPAL Collaboration decided to compare fixed orders 
in perturbation theory to experiment to determine the value of the strong coupling at the Z mass, 
as{Mz), with the results shown. Near the two-jet limit, however, the distributions peak, and 
then show a rapid decrease. This behavior cannot be reproduced by a few orders in perturbation 
theory, but requires an all-orders analysis. 
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Figure 12: Distributions of various event shapes as measured by the OPAL Collaboration. 
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As we've already observed, both the thrust distribution da/dT and its MeUin moments 
jT^a. dTT^ da/dT are infrared safe. The rapid turnover near T — 1 may be traced to a sin- 
gularity in the thrust distribution there, which appears at 0{as) as 



1 daiQ) _ 5(i_j.)_2C^«.(Q) 



ln(l -T) 



+ 



(To dT 71 

IfdTT^^ = l-C,'^ln^N+.... (78) 

(To Jo dT TT 

Here the "plus" distribution [...]+ is associated with the sum of virtual-gluon and real-gluon 
corrections. It will be defined below. 

Prom a physical point of view, taking T — > 1 is like taking energy resolution to zero in QED, 
and we expect the cross section to vanish in that limit, simply because radiation is required when 
charged particles arc accelerated. In QED this is generally more a matter of principle than of 
practical concern, because the coupling is small. The strong coupling even at the Z mass is not 
so small, however, and the corrections of (78) are not only significant, they dominate the cross 
section near the peak. This is characteristic of event shape distributions, for which the bulk 
of events are precisely in region where corrections from every order in perturbation theory are 
large. This sets the stage for a discussion of resummation in QCD, which enables us to bring such 
corrections under control. Wc will come back this issue later, after a discussion of applications 
of perturbative QCD to deep-inelastic scattering. 



4 Deep-inelastic Scattering and Factorization 

In this section, we review the basic structure of deep-inelastic scattering, showing how QCD at 
once justifies and generalizes the parton model. We will introduce the method of factorization, 
which enables us to use infrared safety, and hence asymptotic freedom in DIS and other hard 
scattering processes with hadrons in the initial state. In this section, we will deal primarily with 
low-order examples. A brief discussion of what is involved in proofs of factorization theorems 
will be given in the final section. 

We'll begin with a quick review of DIS in the parton model, followed by a statement of its 
generalization to QCD in the form of factorized structure functions. We go on to sketch the 
procedure of calculating perturbative corrections to the structure functions, including the deriva- 
tion of one-loop evolution kernels (splitting functions). All these calculations are elementary, but 
they serve as a prototype for the analysis from which we derive predictions for both new-physics 
signals and QCD backgrounds and tests. 



4.1 DIS and the parton model 

As introduced in Sec. 1, deep-inelastic scattering (DIS) is the reaction e{k) + N{p) — > e{k') + 
Xhadronic, whcrc thc momentum transfer q — k — k'is large. Fig. 13 represents the basic process. 
The corresponding (unpolarized) cross section can be broken down into the product of leptonic 
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Figure 13: DIS kinematics. 



and hadronic tensors, in just the same way as e"'"e annihilation, 

d^h' 1 

^ ^ W ' ^^^^ 

with the leptonic and hadronic parts given by lowest-order QED and by hadronic matrix elements, 
respectively, 

W^JJ" ^ ^ E j:{N{p,a)\UO)\X) 

spins a X 

x{X I J,(0) I N{p,a)){27i)^d\px -q-p). (80) 

The superscript reminds us that we are neglecting Z boson exchange. The symmetries of the 
strong interactions, already mentioned in connection with Eq. (4) for the DIS cross section, lead 
to the following decomposition of W^" in terms of invariant tensors, 

+ (p, + % {£)) {p. + (^)) W^ix, q') , (81) 
which are related to the dimensionless structure functions of Eq. (4) by 

F^ix, Q') ^ F^ix, Q')=p-q Wr , (82) 

where again, the "scaling" variable is given by x = —q^/2p ■ q = Q^/2p ■ q, and as noted above, 
scaling is the statement that the structure functions depend only on x, not on Q. Contemporary 
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Figure 14: HI data for the structure function F2. 



data in Fig. 14 from the HI experiment at DESY shows that for a wide range of x and Q scahng 
is an experimental, albeit approximate, statement of fact. Recalling the discussion of Sec. 2, the 
parton model cross section for electron-nucleon DIS at a given x and is given by the product 
of the lowest-order electromagnetic "Born" cross section for electron-quark scattering, times the 
probability distribution for quarks (of flavors /) within the hadrons at momentum fraction x, 

cia(^^)(p, g) = E /' da^,J'f\^p, q)cPfMO ■ (83) 
/ ° 

In the parton model, the quark distributions, (pf/N^x) are independent of Q. 

We can break down the parton model cross section (83) into leptonic and hadronic tensors, 
with the latter given by 

From here we identify structure functions in the parton model, following the general analysis 
(80), (81) and (82), which gives 

F^ix) = diF^,f^'\x/i)<Pf/^{0 
f 

Fn^) = j:t^F^'^'\x/o<i>f/Nio, (85) 
/ ° ^ 

where the Born structure functions are found by the replacing the external hadron by external 
quarks. In Eq. (85), the Born structure functions are then 

2F7^(°)(^) = F^^'-'\z) = Qf^6{l - z) . (86) 
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Inserting (86) into (85), we find 

2xFr{x) = F^f^'^ (x) ^Y.Qf^ h/N{x) , (87) 

/ 

which is known as the Callan-Gross relation [43] . The structure functions are thus proportional 
to the quark distributions in the parton model, and hence to each other. As anticipated in Sec. 
2, this relation depends on the quark spin, and is thus analogous in significiance to the 1 + cos^ 6 
dependence of jets in e"'"e~ annihilation, providing direct insight into parton spin. 

4.2 Factorization in DIS 

We now discuss how QCD incorporates the successes of the parton model. The challenge is to 
use asymptotic freedom in quantities like structure functions. Structure functions clearly depend 
upon hadronic structure and are therefore not infrared safe, so we have no hope of computing 
them in perturbation theory. Factorization is the term that describes a general approach to this 
problem. It is applicable in one form or another for any cross section that includes a short- 
distance subprocess. 

The particular form of factorization that we will discuss is adapted to generalize the parton 
model, and is sometimes referred to as "coUinear" factorization. We'll mention some other 
possibilities later, but coUinear factorization is the best understood and most highly developed. 
For the remainder of this section, we'll use the terms factorization and coUinear factorization 
interchangeably. 

A factorized cross section will be a product or convolution of two pieces: an infrared safe part 
(short-distance), which is calculable in perturbative, QCD, and an infrared sensitive part (long 
distance), which although not perturbatively calculable, is measureable and universal among a 
class of factorizable processes. For DIS and related cross sections, the infrared factors are parton 
distributions, or PDFs. The same general method of factorization can also be applied to elastic 
scattering amplitudes at large momentum transfer [8]. 

The form of a factorized cross section is very close to the parton model, but unlike the 
parton model, the short-distance functions are not normalized uniquely by partonic cross sections. 
Beyond the lowest order in perturbation theory, the split between short- and long-distance is not 
unique. In fact, the parton model would be exact if there were a lower limit to the lifetime 
of virtual states in QCD, corresponding to an upper limit in the energy deficits of such states. 
If this were the case, once the momentum transfer Q were large enough, quantum interference 
between QCD and electroweak scattering would disappear. This is not the case, of course, as 
illustrated by Fig. 15, which shows the emission of a single gluon just before a quark absorbs 
a photon. The lifetime of the quark/gluon virtual decreases with an increase in the transverse 
momentum of the gluon, kr, which is limited only by the available energy. On the other hand, 
for very small kx, the diagram develops coUinear singularities. It is these two regimes we need 
to factorize. These issues will be clarified by studying the first nontrivial order. 



39 




Figure 15: Gluon emission at lowest order in a DIS process. 



4.3 Field theory corrections 

To see how factorized cross sections work, we study DIS not of a hadron, but of a quark. What 
does this mean in perturbation theory? To compute a quark structure function, we must work 
in the (dimensionally) regulated theory, with the number of dimensions n > 4. Although, as 
emphasized above, this is not the real theory of QCD, infrared safe quantities computed in this 
theory can be returned to four dimensions, and incorporated into realistic QCD. So let's go 
ahead. 

4.3.1 Factorization for quark DIS 

Let us begin with the statement of factorization, in this case for the structure function of a 
quark, given by the convolution of a short- distance, infrared safe function, C, which generalizes 
the parton model, and an infrared- sensitive function 0, 



The first equality defines the convolution indicated by the second. The sum is over all flavors of 
quark and antiquark, as well as the gluon. Radiation may produce any of these partons in the 
initial state, and the transition from long to short distances may occur for any parton. We note, 
of course, that when i = g, the coefficient function must be at least order as, since the gluon must 
fluctuate back to a quark pair to absorb a photon. The variable fip is the factorization scale, 
which will separate short- and long-distance dynamics, while fi is the renormalization scale. For 
our purposes, it is convenient to identify these two scales, fip = fJ', although we need not do so. 
We also often encounter the direct choice fi = Q, which emphasizes the perturbative calculability 







of C, 
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Eq. (88), however, illustrates the general statement of factorization. 

The zeroth order expansion of the structure function F2 is found from the Born graph of 
elastic quark scattering, while NLO contributions to Eq. (88) are calculated from the diagrams 
of Fig. 16 in n dimensions. We suppress, but hopefully remember, dependence on the variable 
e = 2 — which will appear as poles in the infrared-sensitive factors of this expression. 




Figure 16: NLO diagrams for F2. 
At lowest order (no QCD), the distribution of a parton in a parton is a delta function, 

</,'(0 = V 5(1-0, (90) 

which expresses the obvious fact that in the absence of any interactions, the quark's momentum 
fraction is conserved. Then, comparing the factorized structure functions of Eq. (89) to the Born 
quark structure functions in Eq. (86), we confirm that at zeroth order the factorized coefficient 
functions are the Born structure functions (86) 

clT (1) = Ql - ^/O , (91) 

so that at lowest order the factorized quark structure function reduces to 

= Q) ['d^6{l-x/0S{l-0 

J X 

= Q}xS{l-x), (92) 
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exactly as in the parton model. We are now ready to incorporate QCD corrections. 



4.4 Quark DIS to one loop 

Moving on to NLO (often termed one loop, even though gluon emission at this order is really 
tree level), as shown in Fig. 16, we recognize that we are going to have contributions from two- 
and three-particle phase space. Recalling our discussion of infrared and coUinear divergences, 
we anticipate that there will be at least some cancellation of poles in e. The treatment of these 
technical issues is simplified by the introduction of plus distributions. Like delta functions, plus 
distributions are generalized functions, and are defined by their integrals with smooth functions. 
The definitions of the plus distributions we'll need are, 

Jo (1 — x)+ Jo (1 — x) 

with an obvious to generalization to higher powers of logs and other numerator functions. In 

factorized cross sections, the role of f{x) will be played by products of parton distributions and 
other smooth functions. At first order, the f{x) term corresponds to the emission of a real gluon, 
with momentum fraction 1 — x, while the /(I) term will correspond to virtual gluon corrections, 
with elastic kinematics for the quark. A very special distribution, one we will encounter shortly, 
is is the evolution kernel or splitting function for quark-to-quark evolution. 



1 + x^ 



X 



(94) 



For us Fqq will have the interpretation of the probability of gluon emission per unit log transverse 
momentum. 

Let us now relate the first two orders (LO, NLO) on the left- and right-hand sides of Eq. (88). 
The left-hand side is defined uniquely by the diagrams in (regularized) perturbation theory, but 
we need a criterion for splitting the result from the left into a coefficient function and parton 
distribution on the right. Such a criterion is usually called a factorization scheme. In any scheme, 
we can expand the right-hand side in orders of ag, 

F-^ (x, g^) = C^^ ® + ^ c^^'^) + ^ ® + . . . , (95) 

where we introduce the notation /(a^) = Z]n(<^s/27r)"/'^"). Into this template we can insert the 
explicit calculation of the left-hand side, given for F2 by 

F^'^{x,Q') = Q}xl 6{l-x) 
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and for the correction to the Callan-Gross relation by 

2xFr'{x,Q') 



a. 



+ 



X, Q') -Cf — x' 

TT 



(96) 



(97) 



In Eq. (96), the dots indicate additional £-dependence that vanishes in four dimensions. The 
remaining integral over kr refers to the transverse momentum of a virtual or real gluon in the 
diagrams of Fig. 16. This term illustrates the interpretation of the evolution kernel (94) as a 
probability for gluon emission per unit logarithm of kr-, in four dimensions. At low, kx-, however, 
we see that infrared regulation is required, as confirmed by the necessity of £ < 0, i.e., n > 4 to 
define the integral. 

The calculation that leads to Eq. (96) is straightforward, and a clear exposition of all details, 

for both quarks and gluons may be found in [44]. The two- and three-body final states both 
contain double poles in £, which, however, cancel in the inclusive hadronic cross section. The 
factor in the final line contains the remaining infrared sensitivity, which can be identified as a 
coUinear divergence. The factor is here the renormalization scale in the dimensionally regular- 
ized theory and the factor 47re~'''^, with 7^ the Euler constant, is characteristic of all one-loop 
integrals. As expected, despite the cancellation of double poles, this structure function cannot 
be returned to four dimensions. We are now ready to discuss schemes for choosing the infrared 
safe coefficient functions. 



4.4.1 Factorization Schemes 

As long as £ < 0, we are free to perform the k"^ integral in Eq. (96), 



/7P 1 

4+^^ -e^ ■ ^ 



The pole generated by this integral, times the DGLAP sphtting function, is the only term in Eq. 
(96) that must be absorbed into the definition of 0*^^^ in Eq. (95). In the minimal subtraction, 
or MS, scheme for the quark distribution in n = 2 — 2£ dimensions, this is exactly what we do, 

MS : <t>^^l{x, I,') = {ATTi^'e-^-y P,,{x) £ ^ , (99) 

with the factorization scale. The advantage of this choice is its technical simplicity and gener- 
ality. The coefficient function at one loop is easy to find, 

{x)y^g = Pqq{x) ln(g7/^^) + /^-Independent , (100) 
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where all of its factorization scale dependence is specified by the splitting function. Of course, 
other choices of scheme are possible, by shifting other finite terms from the coefficient functions 
into the distributions. Closer to the parton model than the MS scheme is the DIS scheme, which 
is defined to absorb all corrections into the distribution for deep-inelastic scattering, 

DIS : = Fr^'\x,f,'). (101) 

In this case, 

Ci'\x)ms = xP,,(x)ln(QV/x2) +0. (102) 

Notice that even in DIS scheme, the coefficient inherits the same dependence on the factorization 
scale as in the MS scheme. Other schemes, closer to other processes are possible, but nowadays, 
most calculations are presented in a minimal scheme. 

4.4.2 Using the regulated theory: parton distributions for real hadrons 

As we have repeatedly observed, IR-regulated QCD is not the real QCD, but since it differs 
only in its infrared sensitivity, we should be able to use infrared safe functions like the Ca^'s 
even in four dimensions, where we cannot calculate parton distributions for real hadrons, at least 
perturbatively. This will enable us to get parton distributions from real hadrons. Here is how it 
works: 

• Compute F^'^, F^^ (for weak bosons as well). 

• Define a factorization scheme; find infrared safe C's. 

• Use these coefficients in the full theory for hadron A^, 

= CT®(t>a/N. (103) 

• Measure the physical structure functions Fi (not neglecting W and Z boson exchange and 
neutrino scattering); then use the known C's to derive the full set of 0a/iv for q = q,q, G. 

Compared to our discussion above, multiple flavors and structure functions essentially com- 
plicate the technicalities, but not the logic. Of course, the final step is a highly nontrivial exercise 
in undoing the convolutions that relate distributions to observables. Assuming we can carry out 
this procedure, however, the result is a collection of parton distributions, 0a/iv(C, A*^), for all 
partons a (quarks and antiquarks of various flavors, and the gluon) in the nucleon (proton and 
neutron). The parton distributions can now be used in any process that factorizes into distri- 
butions times pcrturbative hard-scattering functions. What we have not explained yet is how to 
combine different choices of /x^. To that end, we turn to evolution. 



44 



5 Evolution 



Evolution enables us to determine the Q^-dependence of structure functions, and more than that, 
to extrapolate parton distributions from one energy scale to another. It is the cornerstone of 
predictions at high energy both for QCD and new physics signals. 

For the purposes of this discussion we set the factorization scale equal to the renormalization 
scale, hf — A*, so we write Ca {x/^, Q'^f/j,'^, asi/j)). We now choose ji — Q^so that F2, for example, 
is given for hadron A by 

Fl\x, Q') = [ Cr (1 1, as{Q)^ (t>a/A{^i Q') . (104) 
The larger is Q, the better the perturbative expansion, 

cr (|, 1, «s(g)) = E (^)" ^2"^"^ (f ) • (105) 

To make this work for us, we only need the distribution aX. 11 = Q. 

Assuming that we have determined our PDFs at some scale Qo, evolution makes it possible 
to determine them at any other scale Q, and hence to determine the structure functions (for 
example) at any other Q. Our only requirement is that both Q and Qq should be large enough 
that both as{Q) and otsiQo) are small. In particular, extrapolations to arbitrarily large energies 
are possible. 

We can illustrate these ideas most easily by studying a "nonsinglet" structure function, defined 
as the difference between proton and neutron structure functions 

^7NS _ _ ^7n ^ (106) 

which may be factorized in terms of a single nonsinglet parton distribution, 

Ff\x,Q') = Y.tdiCr^(-.'^,as{^)\ct>^s{i,^'). (107) 

In this expression, the contributions of gluons and antiquarks have cancelled, and the nonsinglet 
distribution is a sum over quark flavors. In fact, these cancellations are relevant only beyond one 
loop; at one loop C^^ = Cf. 

It is now convenient to take the Mellin transform of Eq. (107), defined by 



fiN) ^ j dx x^-^ fix) . (108) 
This has the advantage of reducing convolutions of this kind to products, 

f{x) = I'dyg (-\ h{y) ^ f{N) = g{N) h{N + 1) . (109) 
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In our case, we find 



F,^^^(iv,g2) = cf' (n, 0Ns(iv,/x^) , (110) 

where in a slight shift of notation we define 0ns (-^, 1^'^) = lo A*^) here. 

We now observe that F2^^{N,Q^) is physical, so that, as in Eq. (12), it cannot depend on 
the factorization/renormalization scale, /i, 

f^-^F^'''{N,Q') = 0. (Ill) 

Applied to Eq. (110), this means that variations in the coefficient function and in the parton 
distribution must compensate each other, 

11-^ ln0Ns(A^,A*^) = -7ns(A^, Qis( A*)) 

7Ns(iV, «,(/.)) = i^^lnCf^{N,as{i^)) . (112) 

The function 7ns (-^, ctsifJ^)) plays the role of a separation constant. It is a function of ag and N 
only, because these are the only variables that C and hold in common. In particular, it has 
an infrared safe expansion in a^. Thus, to get the 7's we need to know only C's, which we can 
get from the infrared regulated theory as above. The anomalous dimensions 7(A^) found in this 
way, however, enable us to determine the iip dependence of the parton distributions. For any 
factorized cross section, we may choose /ip at the hard scale in the problem, and make predictions 
at high energy based on observations at lower energy. 

In summary, the Q-dependence of factorized cross sections is determined by perturbation 
theory through evolution. This procedure was how we found out that QCD is the "right" theory 
for the strong interactions. It is also the way QCD predicts physics at as-yet unseen scales. 

5.1 Evolution at one loop 

We now have all the pieces necessary to derive the consequences of asymptotic freedom and 
factorization in deep-inelastic scattering. According to (112), we can compute the anomalous 
dimensions directly from the coefficient functions, noting that the answer is the same in either 
(or any) factorization scheme, 

= 11^ { {as/2Ti) Pgg{N) ln(gV/x') + indep. } 



d/i 

TT Jo 
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7ns- 



(113) 



We have used the exphcit form of the nonsinglet sphtting function, Eq. (94). These nonsinglet 
anomalous dimensions are the same as those for the full structure functions, to lowest order. To 
confirm the integrals the identity 



(l-x^)/(l-x) 



k-l 
i=0 



(114) 



may be helpful. 

Having done the work, it's just a small step to some extraordinary physics predictions ^. 
Going back to Eq. (112), the very anomalous dimensions computed from the NLO coefficient 
functions control the factorization scale dependence of the parton distributions. 



(pNsiN, iji^) = -^m{N,as(iJi)) 0ns(-?V,A«^) , 



(115) 



The solution to this equation is an exponential, where the exponent is an integral over the scale 
appearing in the running coupling. 



(pNs{N,i/) = 0Ns(^^,Aio) exp 



1 f/^' 

— 7Ns(A^,tt.(/^)) 



2 ^Mg 



(116) 



Of course, we can evaluate this expression for any theory, once we determine the anomalous 
dimensions and the running of the coupling as above. If (the generic case) the theory is not 
asymptotically free, the coupling increases rather than decreases with /i', on the basis of its 
one-loop beta function. At the same time the anomalous dimensions, the analogs of 7^^g(A^, a^), 
remain positive. Once the coupling is large, of course, it may no longer be possible to use the 
one-loop anomalous dimensions, and it is still possible that the coupling constant might approach 
a finite constant (a fixed point) due to higher-order terms in the beta function. This, however, 
would predict violations of scaling by a power in fi, and hence in Q [17]. In contemporary terms, 
for such theories the parton distributions decrease rapidly in /x, and hence the moments of the 
structure functions decrease very strongly with Q. In physical terms, a decrease in the moments 

^As the written form of these lectures were being prepared, the 2004 Nobel Prize honoring the discovery of 
asymptotic freedom was awarded. The field theoretic explanation of scaling and scale breaking, in much the same 
manner as presented here, was the original link of QCD to experiment and remains one of its strongest links to 
this day 
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corresponds to a rapid shift of the parton distributions from larger to smaller x as Q increases. 
This was the paradox of scahng in 1973: all known field theories seemed to predict too much 
radiation for partons to retain their momentum fractions x at large Q^. 

In QCD, of course, asymptotic freedom changes this. The QCD coupling (11) decreases with 
scale. We can perform the /i' integral in the exponent of (116) easily, to get 

^ ^..(iv,«s)(;^|gg)""", (117) 



which is also often written as 




cl>m{N,Q') = 0Ns(A^,Q^) ■ (118) 



Scale breaking has not disappeared, but it is now logarithmic, and hence very mild once Q is 
significantly larger than the scale Aqcd- This establishes both the consistency of QCD with the 
parton model, and a constellation of predictions for improvements to the parton model, which 
have been confirmed in the intervening years. 

The scaling violation predictions of QCD are easiest to compute via moments of the parton 
distributions as above. For numerical evaluation, however, it is often more convenient to invert 
the moments, which takes the simple rate equation for the parton density in (112) back to a 
convolution 

in terms of the very same splitting function, 

dx x^-^ Pgg{x, as) = --fgg{N, «,) + ... = x^'^ Pm{x, «,) + ..., (120) 

where again, the nonsinglet and quark-quark splitting functions are the same to lowest order. 
Equation (119), commonly referred to as the DGLAP evolution equation [46], also provides a 
nice physical interpretation of evolution, as a sequence of radiation processes. At each such step, 
the parton momentum fraction ^ is degraded to x < ^, in an essentially probabalistic manner. 



5.2 Beyond nonsinglets 

The generalization of these results to the full system of parton distributions is simple, at least in 
principle. As long as the sequential radiation processes that drive evolution involve only gluons, 
the evolutions of individual quark flavors are independent. The (singlet) structure functions for 
individual hadrons, however, receive contributions from fermionic as well as gluonic radiation, 
which couples the evolutions of all the quarks, antiquarks and the gluon. Compared to Eq. (119), 
the only change is that the nonsinglet splitting function is replaced by a matrix P^a, each element 
of which describes the production of partons of flavor b from those of flavor a, 

d /"i d^ 

11-^ (j)b/A{x, fl^) = - / ^ Aa(C,«.(/^)) 0aM(C,/^^) , (121) 
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or in moment space 

^l^^h/A{N,^l^) ^ - 76a(iV,«s(/x))0aM(iV,/^'), (122) 

^ b=q,q,G 

in terms of a matrix of anomalous dimensions, jbaiN). Of special interest is the low- a; behavior 
of Pgg{z) 1/z, whlch corresponds to ^gg{N) ~ 1/{N — 1). For N 1 moments, then, the 
effects of evolution can be large, simply because large anomalous dimensions raise the logarithms 
corresponding to those in (118) above to a high power. The strong scale-breaking effects of such 
terms in the anomalous dimensions are clearly seen in the data of Fig. 17. The systematic study 
of \ow-x evolution is a subfield in its own right, which starts with the famous BFKL equation 
[47], which organizes the leading logarithms in x at each order of ccg. An introduction can be 
found in Ref . [35] , and a summary of recent highlights in Ref . [48] . 

5.3 Summary 

Here is a summary our discussion of factorization and evolution. In the parton model, (t>a/A{x) 
denotes the density of partons a with momentum fraction a distribution that is assumed to be 
quantum-mechanically independent of the hard scattering at momentum transfer Q, and hence 
may be treated as an independent probability. In QCD, 0a/^(a;, /x) represents the same density, 
but only of partons with transverse momentum < /x, as illustrated by the integral of Eq. (96). It 
is only these partons whose production may be considered incoherent with the hard scattering. 

If there were a maximum transverse momentum Qq for partons in the nucleon, 0(a;, Qq) would 
freeze for // > Qq, and the theory would revert to the parton model above that scale. This is never 
the case, however, in a renormalizable field theory, and scale breaking measures the change in 
the density as the maximum transverse momentum increases. Of course, the structure functions 
and cross sections that we compute still depend on our choice of through uncomputed higher 
orders in C and evolution. At present, evolution for DIS is routinely implemented fully up to 
NLO, which requires order a1 in the splitting functions Pab- This formalism can successfully 
describe the behavior of structure functions over a wide range of x, as illustrated by Figs. 14 and 
17. The full sphtting functions up to have just become available within the past year [49], 
and the development of a phenomenology at NNLO is underway [50] . 

6 Applications and Extensions 

In this final section, we will touch on a (very incomplete) subset of the generalizations of the jet 
and DIS analyses we have just described, including factorization for hadron-hadron scattering, the 
basis of all-orders factorization proofs and of resummations, and a few current areas of progress 
relevant to the LHC program. 
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Figure 17: DIS data at low x, compared to fits from [45]. 
6.1 Factorized hadron hadron cross sections 

Space does not allow an extensive description of the application of factorization and evolution 
to hadron hadron scattering [30, 35], but even a cursory description shows the utility of these 
techniques. The factorization formalism can be extended to the semi-inclusive production of 
heavy particles or systems of jets, labeled F, of mass Q. Indeed, the decay of a system of heavy 
particles, be it a top quark pair or SUSY signal, usually results in hadronic jets. An essential 
element of new-particle searches is to tease the signals of these decays out of the QCD and 
electroweak background. 

The production of final-state F requires the short-distance collisions of a parton from each 
hadron, and the same pinch surface and power counting analysis as for e^e~ annihilation and 
DIS shows that the entire leading power in Q comes from only a single (physical) parton from 
each hadron [30]. The essential requirement for factorization is once again that we sum over 
states involving the emission of soft particles and/or the rearrangement of collinear radiation. 
It is important that the definition of the semi-inclusive final state not involve any scales much 
smaller than Q, the mass of the observed final state, because the largest corrections will generally 
be an inverse power of the smallest such scale. 

A factorized cross section for hadron-hadron scattering describes the production of final states 
by a direct generalization of DIS factorization, as a convolution of two parton distributions with 
a hard-scattering function, 

(JAB-.F{Q){Q) = Yj J dxdy4>a/A{x,^) (f)b/B{y,fJ')^ab^FiQ){xpA,yPB,Q,fJ',(ys{fJ')) 
a,b=q,q,G ^ 
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(123) 



where we have exhibited a theta function to emphasize the requirement that the partonic system 
has sufficient energy to produce system F of mass Q. As usual, the hard-scattering function 
a begins at lowest order with a Born cross section, and depends on the choice of /j, through 
uncomputed higher orders. 

An illustrative 0{as) (NLO) hard-scattering cross section is the quark- antiquark annihilation 
cross section aqq for the Drell-Yan process (qq — > 7 or Z — > given in MS by 



4? = <TBorn(g^) |p,,(^)ln(^^j+2(l + z^) 



ln(l - z) 
1-z 



(l-z) 




(124) 



where z = Q'^/xaXi,S with a/S" the total center of mass energy squared, and where o"Born describes 
the electroweak annihilation of q q. We see again a plus distribution, whose presence reflects the 
cancellation of real and virtual soft gluons. 

A perturbative prediction for this process, based on factorization, is compared to data at high 
energy (7 and Z), is shown in Fig. 18 from CDF at the Tevatron, The peak due to the Z is clearly 
visible, and indeed this was the way the Z was ffist directly seen. Also shown is data for the 
forward-backward asymmetry associated with the parity-violating predictions of the standard 
model, and some extenstions of it. Notice that this asymmetry, with or without a new physics 
signal, is calculable precisely because it is part of the hard scattering function. Such tests, for 
signals of new states or violations of symmetries, are the prototypes for the detection of new 
physics in hadronic scattering. 

The successes and remaining uncertainties of predictions based on factorizcd perturbation 
theory are well illustrated by the data for inclusive high-]?^ jet cross sections, based on formulas 
that are analogous to those for Drell-Yan, but (much, much) more elaborate even at NLO [52]. 
The data track perturbative predictions over many orders of magnitude, but even between the 
two Fermilab experiments, differences in analysis of the still-preliminary data [53] show that there 
is still work to be done. 

The LHC program for the discovery of new particles and/or new strong interactions at high 
energy relies heavily on the factorization formalism. So long as it involves heavy states, each 
prediction of an extension of the standard model appears in the hard scattering function of the 
factorized cross section. The parton distributions are the same as those observed in DIS at lower 
energy, now extrapolated by evolution. 



6.2 Factorization proofs and matrix elements 

In the previous section we have made extensive use of factorization for the DIS structure functions, 
and now we have seen how they generalize to hadronic collisions. Let us come back briefly to the 
proof of the forms (103) for the former and (123) for the latter. 
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Figure 18: Drell-Yan data from the Tevatron [51]. 



For deep-inelastic scattering, there are a number of ways to provide a proof. One way is 
based on the hght-cone expansion of Eq. (14), and is most conveniently formulated in terms of 
moments. 

Alternatively, we may begin wiith the pinch surfaces and hence physical pictures for the 
process eA^ — e + X, where we assume a large momentum transfer, Q. These are shown in Fig. 
21, and generally involve jets in the final state, in addition to collinear fragments of the incoming 
nucleon. The figure shows a single parton, labeled i combining with the off-shell photon of 
labeled q to produce a final state. Physical pictures with more collinear partons are possible. As 
in the discussion of jet production in e"'"e~ annihilation above, at leading power in Q, they can be 
only unphysically-polarized gluons (in covariant gauges). It is worth noting that the requirement 
of a physical picture implies that the momentum fraction of the parton that initiates the hard 
scattering obeys 1 > ^ > x, for much the same reasons as in our low-order examples above; all 
the partons in a jet or hadron must be moving in the same direction. For the inclusive DIS cross 
section, we sum over all final states, and information on the time development of the final state 
jets is lost, just as in the total e"'"e~ annihilation cross section. For DIS, however, the initial-state 
evolution that provides the "active" parton (the one that initiates the hard scattering) remains. 

The DIS argument regarding the sum over states has the same content as the optical theorem 
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Figure 19: Recent preliminary CDF jet data [53]. The data is presented here as the ratio of the 
difference between data and the NLO prediction, normahzed to the prediction. As the following 
figure shows, the absolute value of the data is changing by many orders of magnitude over this 
range. 



for DIS: the total 7* p cross section is related to the 7*p forward scattering amplitude. In terms 
of the hadronic tensor of Eq. (81), this can be written as 

= 2ImT^,, (125) 

where the tensor T^y is the expectation of the time-ordered product of currents, 

= ^ / ^'^e^'" < ^iP)\ T Mx)Jy{0) \N{p) > . (126) 

Because of the optical theorem, the long-distance behavior of the hadronic tensor is specified by 
the pinch surfaces of this matrix element, which can be described as the amplitude for forward 
Compton scattering of a nucleon by an off-shell photon. The pinch surfaces in this case are much 
simpler than those for individual final states. The requirement of a physical picture shows that 
no jets are possible aside from the "forward" jet of fragment partons from the incoming nucleon, 
which reform after the hard scattering into an outgoing nucleon of the same total momentum^. 
The hard scattering is reduced to a point, because, as in the e"'"e~ case, on-shell particles that 
travel a finite distance from the hard scattering can never recombine. At the pinch surface, the 
hard-scattering function H depends only on the non- vanishing longitudinal component of the 
active parton's momentum. Once we have reached this stage, the remainder of the factorization 
argument is a systematic disentangling of the upper and lower parts of Fig. 22, in terms of color, 
spin and momentum. A hierarchy of subtractions [30] effects this separation. 

related physical off-diagonal amplitude [9], 7* + p 7 + p called deeply-virtual Compton scattering, offers 
complementary information on hadronic structure, and may also be factorized [54]. 
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Figure 20: Recent preliminary DO jet data [53]. 



The parton distributions factorized in this mannner take the form of the expectation values 
of operators relatively on the light cone. To specify these matrix elements, it is convenient to 
choose the momentum of incoming nucleon h in the plus light-cone direction, and to define an 



opposite-moving light-like vector. 



E 



dy 
2-K 



-i^p+y 



The parton distribution is then 



< h{p, a) I g(0+, y-, 0^) -n ■ 7 g(0) | h{p, a) > , (127) 



where the average over spins for an unpolarized cross section is exhibited. This matrix element 
can be defined in n ■ A = = A'^ gauge. Choosing the incoming hadron as a quark, and 
computing (127) at order a^, we find precisely the MS quark distribution discussed in Sec. 4 
above. The formalism can readily be extended to covariant gauges, and the matrix element 
made gauge invariant, by connecting the quark operators with ordered exponentials (also known 
as nonabelian phases, gauge links and Wilson lines). 



q{y ) n ■ 7 g(0) q{y )Pexp 



ig I din- A{ln^') 



n 



7 g(0) . 



(128) 



These and related forms of parton distributions, including their generalizations to transverse 
momentum distributions, were discussed extensively in Ref. [55]. 

An extensive discussion of proofs of perturbative factorization for semi-inclusive hard scat- 
tering cross sections in hadron hadron scattering [56] may be found in [30], with additional 
comments in Ref. [34]. Here, I will only recall a few of the heuristic arguments. 

The essential challenge when there are two hadrons in the initial state is to show the univer- 
sality of the parton distribution functions between DIS and hadron-hadron scattering. Consider 
a Drell-Yan cross section, in which a quark from one hadron annihilates with an antiquark from 
the other to form a virtual photon, W or Z (generic mass Q). We'll think of a range in kinematics 
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Figure 21: Pinch surface for DIS. 

where the quark and antiquark carry a large fraction x of their parent hadron's momentum, so 
that there is no confusion about which of the pair comes from which hadron. 

It is natural to wonder whether, with two hadrons in the initial state the color fields of 
each hadron might permeate the other prior to the annihilation, modifying the colors, or even 
transverse momenta, of the annihilating pair in such a way as to increase or decrease their 
amplitude for the annihilation. A proof of factorization must show that such effects decrease 
as a power of the pair invariant mass, Q. The technical arguments rely on the formalism we 
have developed above in Sec. 3. The summation over final states that defines the semi-inclusive 
cross section results in an expression for which the only pinch surfaces that remain logarithmic in 
their power counting after the sum over final states are in one-to-one correspondence with those 
found in the inclusive DIS cross section, separately for the two hadrons. The hadron-hadron 
cross section is then factorizable in the same manner as deep-inelastic scattering. 

The physical reason for this rather simple result is not far to seek. The Lorentz contracted 
fields of incident particles do not overlap until the moment of the scattering [57]. As a result, 
initial-state interactions that couple the two hadrons disappear at high enough energies. All 
remaining initial-state interactions are internal to the hadrons, and are specified by the same 
parton densities as in deep-inelastic scattering. It is a matter of relativistic causality: as the 
relative velocities of the two hadrons approaches the speed of light, they are unable to exchange 
signals of any kind. Correspondingly, the annihilation process itself (or other hard-scattering 
reaction) must also occur on a very short time scale, and final-state interactions, after the hard 
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Figure 22: Pinch surface for forward Compton scattering. 

scattering, are simply too late to affect large momentum transfer and/or creation of the heavy 
particle(s). 

6.3 Resummations: a closer look at the final state 

These lectures have concentrated on applications of asymptotic freedom, infrared safety and 
factorization in semi-inclusive cross sections with only a single hard scale. A closer look at 
the final state, in terms of shape functions, generally introduces another scale, as noted in Sec. 
3.5. The second scale may be pcrturbative, but ratios of the two perturbativc scales may be 
large at any fixed order in perturbation theory. For example, in the limit of low jet masses in 
Eq. (74) (thrust goes to 1) or of low angularities in Eq. (77), these corrections build up as 
(l/MjJ al' In^""^ (Mji/g) and (l/ea)a^ ln^"~^ e^, respectively. Note the overall additional power 
singularity at zero jet mass in each case. 

In cases like jet masses and angularities, we can control these large corrections to all orders in 
perturbation theory, a process known as resummation. Resummations, like evolution in DIS, can 
usually be derived from a factorization formula, as in Eq. (112) for evolution. We can illustrate 
this process for the angularities (77). Details are given in Ref. [40], but the basic observations 
are simple enough. In the limit of small Ca, the final state consists of two well-coUimated jets, 
accompanied by very soft wide-angle radiation. The soft radiation cannot resolve the internal 
structure of the jets, and is equivalent to the soft radiation emitted by a pair of point-like 
recoil-less sources in the fundamental representation of the gauge group (quark and antiquark 
representations). In this limit, the cross section factorizes as shown in Fig. 23. Here each of the 
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Figure 23: Factorization near the two-jet limit. 



individual factors can be given a field-theoretic interpretation, an example being the functions 
associated with the jets, 
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(129) 



As in the case of the matrix element that defines the quark distribution (127) with (128), we 
have introduced an ordered exponential in the direction of the vector ^ to make the jet function 
gauge-invariant. We also need to introduce a factorization scale, fi, to separate the hard, jet 
and soft functions. Both the scale fi and the vector C,'^ are artifacts of the factorization; neither 
appears in the cross section itself, but the factorization is valid up to positive powers of Cq, which 
is small. The factorized cross section therefore obeys two equations, analogous to (12), 



da 
dfi 
da 
1^ 








(130) 



The solutions to these equations are a generalization of the evolution of parton distributions, 
and as in that case, they are most simply presented in a transform space, in this case a Laplace 
transform with respect to the angularities. The limit of small corresponds directly to large 
values of the transform variable u, dependence on which exponentiates. 
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(131) 



57 



where the exponent -E'(z^) is given by 
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(132) 



In the exponent, the presence of two integrals reflects the two consistency equations (130), and 
the functions Aq, Bq are infrared-safe anomalous dimensions. In particular A = (a^/Tr) Cp + ■ ■ ■ 
also appears in Pqq{x,as) = A{as)/{1 — x) + . . ., and is now known to three loops [49]. For 
quark jets, B = (as/7r)(3Ci7/2) + ... is somewhat more process dependent. Expressions like this, 
with a double integral over the running coupling, were flrst derived in studies of the transverse 
momentum distribution for the Drell-Yan cross section, especially as developed by Collins and 
Soper [58] and Sen [59]. 

The inversion of the transform leads to distributions of the sort shown in Fig. 24 for the mass 
of the "heavy" jet in e"'"e~ annihilation, and for the distribution of Z transverse momenta at the 
Tevatron in Fig. 25. 




Figure 24: Heavy jet distribution. From [60]. 



In both flgures, the cross section has a peak very near the kinematic conflguration of elastic 
scattering: for e+e~, Ca = 0, corresponding to a quark-pair final state, and Qt = 0, corresponding 
to the production of a Z in Born approximation. In both cases, the effect of gluonic radiation is 
to turn a delta function distribution into one that rises to a peak near the elastic limit, and then 
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Figure 25: Transverse momentum distribution for the Z at the Tevatron. From [61]. 

slowly falls off away from that limit. The tail of the distribution away from the peak is given by 
low-order perturbation theory, and the position of the peak can be calculated perturbatively as 
well. Nevertheless, the bulk of the inclusive cross section is concentrated near the peak, and in 
this region, nonperturbative corrections are suppressed only by powers off l/{eaQ)i rather than 
1/Q [37]. This effect can be inferred directly from the Eq. (132) [33]. A sketch of the transition 
from resummed perturbation theory to nonperturbative power corrections can be found in Ref. 
[62]. 

6.4 Higher orders 

Predictions based on factorized cross sections in perturbative QCD require both hard-scattering 
functions and parton distributions. It is an understatement that making progress on both of these 
topics is nontrivial in almost every case of practical interest. The past few years, however, have 
seen important progress, through a growing mastery over higher order corrections, and through 
the improved determination of parton distributions from a growing set of data. Certainly, the 
coming of the LHC has helped drive much of this work, but their own internal logic, and the 
special challenges these projects pose, have surely provided strong motivations as well. 

6.4.1 Toward a two-loop phenomenology 

In effect, the transition from leading order to next-to-leading order accuracy at the level of jet 
cross sections took about a decade, and a similar time scale will have been necessary to realize. 
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for example, full NNLO jet cross sections in time for the LHC Given the experimental value 
of as{mz) ~ 1.2 [63], hard-scattering cross sections at the order of a^(Mz) have a nominal 
accuracy in the per cent range. For some semi-inclusive cross sections at momentum transfers 
in the hundreds of GeV, this is also the scale of many nonperturbativc corrections, including 
those associated with uncertainties in the parton distribution functions (see below). Percent 
accuracy is probably attainable for benchmark cross sections, as a foundation from which to 
understand the complex final states associated with the decay of new, heavy particles, and the 
QCD backgrounds to them. 

The development of a NNLO phenomenology is following the pattern of NLO. For some 
time, fully inclusive NNLO cross sections for e^e~ annihilation, DIS, and Drell-Yan annihilation 
(including Higgs production) have been available [64] . More complex is the combination of final 
states in e"'"e~ jet cross sections, and even more so the factorization necessary for jet cross sections 
in hadron-hadron scattering. Even at NLO the state of the art is 4-jet cross sections for e"''e~, 
and up to 3 jets at hadron colliders [52]. 

Within the last few years, the elements necessary for factorized NNLO cross sections are 
beginning to be assembled in earnest. Progress has been made in NNLO e"'"e~ jet cross sections 
[65], while the amplitude at 0{al) has become available for three jets. The past few years 
have also seen the first NNLO two- loop S- matrix elements for the 2 — > 2 scattering [66], an 
essential step toward jet cross sections in hadron-hadron scattering. Important progress is also 
ongoing toward organizing the singular integrals over phase space [67]. Equally significant is the 
calculation of NNLO splitting functions, the fruit of a decade-long program [49] . 

Many of these acheivements have required developing new methods in the organization of 
large numbers of Feynman diagrams [68] often involving astronomical numbers of terms. Hope 
springs eternal that new mathematical insights [69] may simplify the road ahead, at least for the 
calculation of on-shell amplitudes. 

6.4.2 Parton distributions: global fits and uncertainties 

The determination of parton distributions is a bootstrap process. As described in Sec. 3, they 
must be abstracted from the comparison of sets of data to factorized cross sections, with hard- 
scattering functions computed up to a given order. For the past ten years, the calculation has 
been at the level of NLO. Any set of distributions determined this way is limited not only by the 
accuracy of the calculations, but also by the sensitivity of the data to given distributions. For 
example, in DIS, photons or weak vector bosons couple to the quarks at zeroth order, while the 
gluons get into the act only at order in the hard scattering. Naturally, we might expect DIS 
to be less informative about gluons than about quarks, and we look for other processes, like jet 
production, which can be more directly sensitive to the gluon distribution. Other cases where 
new data have provided new and sometimes surprising insights concern the ratios of d-quark and 
u-quark distributions from asymmetry measurements at the Tevatron collider [51], and of 
the d and u from fixed-target Drell-Yan data [70]. 

In any case, once a model set of PDFs has been produced, it provides predictions for other 
factorizable cross sections. As new data become available, it may require us to modify existing 
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models to incorporate new features. This feedback procedure has evolved into the approach of 
"global" fits, which try to use the best available data from a variety of cross sections. (Naturally, 
these are cross sections for which it is anticipated that no new physics signals are confusing 
things!) 

The leading global fits have been provided by the CTEQ [45] and MRST [71] collaborations. 
Sets of distributions are given in terms of the parameters specified in recent work for each partons 
by 

xf{x,Qo) = AoX^i(l-x)^2e^«^(l + e^*x)^5 (133) 
for each of the CTEQ partons, and 

xq{x,Qo) = A{l-xy{l + ex^-^ + -fx)x^ 

xg{x,Qo) = Ag{l-xy'{l + egX^-^ + -fgx)x^' - A4l-xy-x-^- (134) 

for the MRST quarks and the gluon, respectively. Figure 26 shows a broad agreement between 
recent sets, but also the difference in some details. The parameters in these functions are fitted to 




Figure 26: Recent parton distribution fits. 

the data, but have little physical significance in themselves, aside from certain built-in structure, 
such as the vanishing of the distributions at a; = 1. An important development of the past few 
years is the drive to quantify uncertainties within the sets [45], by analyzing the response of the 
fits as the parameters are varied. It is now routine to see uncertainties based on these analyses 
quoted by experiments alongside their more familiar systematic and statistical errors. 

7 Closing Comments 

In these lectures, we have only scratched the surface of the many applications of perturbative 
QCD, let alone the many and varied explorations of its nonperturbative structure. From the 
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former, evolution at small x, and its relation to nuclear collisions, the analysis of heavy quark 
decays, polarized scattering effective field theory approaches, and developments in multi-loop 
calculations are just a few very active areas. I hope, however, that these lectures may provide a 
useful introductory perspective to some of the ideas underlying this vast area of research. 
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